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ABSTRACT 


This  report  is  concerned  with  the  possibility  of  using  undersea 
cables  to  communicate  with  undersea  receivers  at  frequencies  in  the 
ultra  low  and  extremely  low  ranges.  Because  of  the  low  data  trans¬ 
mission  rate  at  these  frequencies,  communication  is  here  understood  to 
mean  the  transfer  of  short  messages  of  high  information  content. 

We  start  by  deriving  theoretical  expressions  for  the  electric  and 
magnetic  fields  generated  in  a  conducting  medium  of  infinite  extent  by 
linear  current  sources  (i.t.,  straight  insulated  current-carrying 
cables)  of  finite  and  semi -infinite  length  and  of  infinite  length  with 
gaps  of  either  finite  or  infinitesimal  size.  Next,  using  numerical 
integration  techniques  and,  in  some  cases,  a  parametric  representation 
to  make  our  results  frequency  independent,  we  compute  representative 
numerical  values  of  the  field  amplitudes  for  the  above  discontinuous 
sources.  For  comparison  with  these  data,  we  compute  an  extensive  range 
of  field  amplitudes  for  the  more  frequently  encountered  continuous 
linear  source  cf  infinite  length  (i.e.,  a  very-long  straight  insulated 
current-carrying  cable);  we  also  introduce  a  parametric  representation 
for  the  fields  produced  by  this  reference  source  and  we  show  that  for 
any  fixed  perpendicular  distance  from  the  source  there  is  an  optimum 
frequency  for  electric  field  generation. 

We  show  that  the  discontinuous  sources  generally  produce  greatly 
enhanced  electric  fields  in  the  vicinity  of  their  open  ends,  as  com¬ 
pared  with  the  continuous  source.  However,  at  large  distances  none  of 
the  discontinuous  sources  we  have  considered  produce  larger  fields 
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than  the  continuous  source,  and  in  most  cases  the  discontinuous  sources 
produce  smaller  fields.  These  results  suggest  that  the  optimum  use  of 
the  discontinuous  sources  (short  lengths  of  current-carrying  cable; 
point  electrodes)  is  for  short-range  undersea  communication,  whereas 
long  continuous  current-carrying  cables  appear  to  have  the  most  promise 
for  long-range  communication  via  their  electric  and  magnetic  fields. 

Finally,  because  even  long  cables  have  limited  communication 

ranges  at  all  but  the  lowest  ultra  low  frequencies,  we  consider  the  use 

of  arrays  of  parallel  long  cables  to  produce  measurable  electric  and 

magnetic  fields  throughout  specific  regions  of  the  sea.  Our  data 

indicate,  for  example,  that  a  series  connected  array  of  10  long  cables 

(each  of  length  L)  laid  on  the  floor  of  a  sea  1  km  deep  and  carrying  a 

1000  A,  1  Hz  current  could  produce  measurable  magnetic  fields  throughout 
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the  sea  over  a  horizontal  area  of  about  75L  km  .  Such  arrays  could  find 
use  in  crucial  areas  where  reliable  undersea  communication  is  of  high 
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Note : 


In  this  report  we  use  the  abbreviation  ULF  (ultra-low- 
frequencies)  for  frequencies  less  than  5  Hz,  and  we  use 
ELF  (extremely -low- frequencies)  to  designate  frequencies 
in  the  range  5  Hz  to  3  kHz. 
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I .  INTRODUCTION 


Shortly  after  the  end  of  the  First  World  War  a  pair  of  articles 
appeared  in  the  scientific  literature  describing  experiments  and  theo¬ 
retical  work  on  the  electromagnetic  fields  produced  in,  on,  and  above 
the  sea  by  submerged  cables  carrying  alternating  current  [Drysdale, 

1924;  Butterworth,  1924],  As  described  in  the  article  by  Drysdale 
[1924],  and  in  greater  detail  in  a  later  article  by  Wright  [1953],  this 
work  was  undertaken  as  part  of  a  British  Navy  project  called  "Leader 
Gear",  which  was  concerned  with  the  use  of  the  cable-generated  electro¬ 
magnetic  fields  for  navigation.  No  further  articles  appeared  on  the 
fields  produced  by  submerged  cables  until  after  the  end  of  the  Second 
World  War,  when  Von  Aulock  prepared  two  lengthy  unpublished  reports  on 
the  subject  [Von  Aulock,  1948,  1953].  With  the  exception  of  a  useful 
summary  of  Von  Aulock's  results  by  Kraichman  [1976]  and  several  particu¬ 
larly  pertinent  articles  by  Wait  [1952,  1959,  1960],  little  research 
has  since  been  carried  out  on  the  electromagnetic  fields  produced  by 
submerged  cables.  This  is  unfortunate,  because  it  is  our  belief  that 
the  fields  could  have  important  appl ications  in  undersea  communication. 
The  primary  objective  of  this  report  is  to  take  the  work  just  described 
through  a  further  stage  of  development  and  thus  provide  an  improved 
theoretical  basis  for  studies  of  the  feasibility  of  the  use  of  the 
fields  from  submerged  cables  for  undersea  communication. 

To  achieve  this  objective,  we  start  with  the  comparatively  well- 
known  expressions  for  the-  electromagnetic  fields  produced  in  a  conduct¬ 
ing  medium  of  infinite  extent  by  a  current  element  and  then,  always 
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considering  a  conducting  medium  of  infinite  extent,  to  systematically 
derive  the  other  known  expressions  for  (1)  the  electric  field  produced 
by  a  linear  current  source  of  finite  extent  (i.e.,  by  a  short  current- 
carrying  insulated  wire  or  cable)  and  (2)  the  electric  and  magnetic 
fields  produced  by  a  linear  current  source  of  infinite  length  (i.e.,  by 
a  long  current-carrying  insulated  wire  or  cable).  In  passing  we  derive 
(3)  an  expression  for  the  magnetic  field  produced  by  the  linear  current 
source  of  finite  extent,  which  is  a  new  result,  and  we  solve  this  ex¬ 
pression  numerically  to  provide  illustrative  magnetic  field  data.  We 
then  expand  from  this  beginning  by  deriving  (4)  expressions  for  the 
electric  and  magnetic  fields  produced  by  a  linear  current  source  of 
semi -i nf ini te  extent,  with  particular  emphasis  on  the  fields  produced 
around  the  end  of  the  source  where  the  current  enters  the  conducting 

medium.  These  latter  results  are  then  used  to  give  (5)  expressions  for 

the  electric  and  magnetic  fields  produced  by  two  aligned,  separated 
semi-infinite  linear  current  sources,  i.e.,  the  fields  produced  in  the 
vicinity  of  two  point  electrodes  immersed  in  the  conducting  medium. 

Finally,  we  allow  the  spacing  between  the  ends  of  linear  sources  in  (5) 

to  become  infinitesimal  and  thus  obtain  (6)  expressions  for  the  electric 
and  magnetic  fields  produced  by  a  linear  current  source  of  infinite  ex¬ 
tent  containing  an  infinitesimal  gap.  As  we  will  show,  these  latter 
expressir:'S  can  alsc  be  derived  by  subtracting  the  electric  and  magne¬ 
tic  field  expressions  for  a  current  element  from  the  equivalent  expres¬ 
sions  for  a  linear  current  source  of  infinite  length. 

In  cases  (4)  to  \6)  above  we  present  representative  numerical 


values  for  the  electric  and  magnetic  fields  that  are  produced.  We  be¬ 
lieve  both  these  numerical  values  and  the  theoretical  expressions  for 
the  electromagnetic  fields  are  new.  We  end  with  a  brief  discussion  of 
the  possible  application  of  the  electric  and  magnetic  fields  to  under¬ 
sea  communication. 

There  is  one  technical  feature  of  our  theoretical  approach  that 
has  general  implications  and  which  is  best  discussed  at  this  introduc¬ 
tory  stage.  It  will  be  noted  that  we  always  assume  an  infinitely  thin 
cross-section  for  the  various  (cylindrical)  linear  current  sources  con¬ 
sidered  in  this  work,  whereas  it  is  obvious  that  these  sources  will 
have  a  finite  cross-section  in  practice.  Technically,  this  assumption 
takes  the  form  of  a  replacement  of  a  volume  integral  in  the  Hertz  vec¬ 
tor  expression  by  an  approximate  line  integral.  The  derivation  of  the 
approximate  line  integral  is  given  in  Appendix  A  and  it  is  shown  that 
the  approximation  is  valid  provided  the  shortest  distance  from  the 
observation  point  to  the  surface  of  the  current  source  is  much  greater 
than  the  radius  of  the  source.  Provided  this  condition  is  met,  our 
expressions  should  give  accurate  values  for  the  electric  and  magnetic 
fields  produced  by  the  linear  current  sources  of  finite  cross-section. 

The  most  important  quantity  to  keep  in  mind  in  this  work  is  the 
skin  depth  <5  defined  by 

6  =  U/uiuc)4  (I  .1) 

for  a  magnetic  (p  ^  pQ)  conducting  medium,  where  we  use  u>  for  the  angu¬ 
lar  frequency,  p  for  the  permeability  of  the  medium  (pQ  is  the  permea- 
bility  of  free  space),  and  a  for  its  conductivity.  At  the  frequencies 
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considered  in  this  report,  electromagnetic  fields  propagate  in  a  con¬ 
ducting  medium  with  a  wavelength  given  by  X  =  2n  6  ,  and  evidence  of  this 
wavelength  will  be  seen  in  most  of  the  figures  showing  numerical  data. 

In  all  cases  these  numerical  data  are  computed  for  sea  water,  where  we 
assume  o  =  4.0  S/m  and  u  =  u0  =  4tt  x  IQ'7  H/m.  Table  1.1  lists  numeri¬ 
cal  values  of  skin  depth  for  these  assumed  sea  water  parameters  and  for 
frequencies  in  the  range  0.001  -  1000  Hz. 

The  theory  developed  in  this  report  assumes  that  the  contribution 
from  the  displacement  current  term  in  Maxwell's  equations  can  be  neg¬ 
lected  in  comparison  with  the  contribution  from  the  conduction  current 
term.  This  is  only  possible  for  frequencies  satisfying  the  condition 
o/a)£>>l  ,  where  e  is  the  permittivity  of  the  medium.  In  the  case  of 
sea  water,  which  is  the  conducting  medium  of  primary  interest  in  this 
work,  the  above  condition  implies  that  the  frequencies  must  be  very 
much  less  than  890  MHz  for  our  field  expressions  to  be  valid.  This 
means  that  the  frequencies  must  be  below  the  microwave  range  (the  low¬ 
est  microwave  frequency  is  usually  taken  to  be  300  MHz)  or,  perhaps 
more  strictly,  they  must  be  less  than  100  MHz.  , 

For  frequencies  satisfying  the  condition  o/u.c>>1  (l.e.,  for  fre¬ 
quencies  less  than  100  MHz  in  sea  water)  it  is  easily  shown  th:t  elec¬ 
tric  fields  propagating  in  a  conducting  medium  are  attenuated  at  the 
rate  of  55  db/wavelength.  This  is  a  very  substantial  rate  of  attenua¬ 
tion  and  it  is  obvious  as  a  result  that  large  ranges  for  communication 
through  sea  water  by  means  of  electromagnetic  signals  can  only  be 
achieved  by  using  low  frequencies.  It  is  for  this  reason  that  we  have 
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Table  1.1  Representative  Skin  Depths  for  Sea  Water 
(os  ■  4.0  S/m) 


Frequency 

(Hz) 

Skin  Depth 
(m) 

Frequency 

(Hz) 

Skin  Depth 
(m) 

1000 

8.0 

1 

251 .6 

800 

8.9 

0.8 

281.3 

600 

10.3 

0.6 

324.9 

600 

12.6 

0.4 

397.9 

200 

17.8 

0.2 

562.7 

100 

25.2 

C.l 

795.8 

80 

28.1 

0.08 

889.7 

60 

32.5 

0.06 

1027.3 

40 

39.8 

0.04 

1258.2 

20 

56.3 

0.02 

1779.4 

10 

79.6 

0.01 

2516.5 

8 

89.0 

0.008 

2813.5 

6 

102.7 

0.006 

3248.7 

4 

125.8 

0.004 

3978.9 

2 

177.9 

0.002 

5627.0 

1 

251 .6 

0.001 

7957.7 
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restricted  our  computations  of  numerical  data  to  frequencies  In  the 
extremely-low  (ELF;  frequencies  In  the  range  5  Hz  -  3  kHz)  and  ultra- 
low  (ULF;  frequencies  less  than  5  Hz)  frequency  bands.  At  these  low 
frequencies  the  wavelength  In  sea  water  becomes  usefully  large  (1.58  km 
at  1  Hz)  and  reasonably  large  communication  ranges  become  possible. 
However,  even  at  these  low  frequencies  the  same  basic  attenuation  rate 
of  55  db/wavel ength  still  applies  and  evidence  for  It  will  be  found  In 
all  cf  our  numerical  data. 


II.  CURRENT  ELEMENT  SOURCE 


Suppose  we  have  a  homogeneous  Isotropic  conducting  medium  of  con¬ 
ductivity  o,  dielectric  constant  c,  and  permeability  y.  Assuming  that 
all  the  fields  vary  with  time  t  as  exp(lut).  Maxwell's  equations  can 


be  written  in  the  form 

dlv  E  =  0  (II. 1) 
div  H  =  0  (II. 2) 
curl  E  =  - i to u H  (II. 3) 
curl  H  =  (a  +  i  u)  e)£  +  J  (II. 4) 


where  J  is  the  Impressed  current  density  and  where  E  and  H  are  the  elec¬ 
tric  and  magnetic  field  intensities,  respectively.  We  wish  to  obtain  E 
and  H  for  a  current  element  Id£.  Immersed  in  this  conducting  medium, 
which  we  will  assume  is  of  infinite  extent. 

In  our  derivation,  and  throughout  this  report,  we  will  make  exten¬ 
sive  use  of  the  Hertz  vector  II  [Panofsky  and  Phillips,  1962;  Kraichman, 
1976],  from  which  the  E  and  H  fields  can  be  obtained  by  using  the  follow¬ 
ing  equations  in  conjunction  with  Maxwell's  equations  II. 1  -  II. 4: 

E  =  -y2n  +  grad  div  n,  (II. 5) 

2 

H  =  -7-*—  curl  n.  (II. 6) 

1  u)p 


Here,  y,  the  propagation  constant,  is  defined  by 


y  =  |i  (d  y  (  o  +  iooc)^ 


(II. 7) 


7 


Substitution  of  Equations  (II. 5)  and  (II. 6)  Into  either  (II. 3)  or 
(II. 4)  yields  the  differential  equation 


curl  curl  n  -  grad  d  1  v  n  +  y  n  = 


(II. 8) 


which  can  also  be  written  as 


y2)n 


-J 

a+i  u)  e 


(II. 9) 


This  equation  Is  called  the  inhomogeneous  wave  equation  and  for  our  un¬ 
bounded  conducting  medium  it  has  a  solution  of  the  form  (e.g.,  Walt, 
1959) 


g  •  (u>°> 

V 

For  a  current  element  source  oriented  In  the  z-dlrection  In  a 
cylindrical  coordinate  system,  as  shown  in  Figure  II. 1,  one  can  easily 
verify  that  the  Hertz  vector  at  the  point  P  has  only  a  z-component  given 
by 


nz(p,z) 


Ida 

4tt(cf+1  oocT 


(11.11) 


where 


8 


Figure  II. 1.  Coordinate  system  and  geometry  used  in  the  derivation  of 
the  electric  and  magnetic  fields  produced  in  a  conducting  medium  of 
infinite  extent  by  a  current  element  Idt. 


Using  this  component  of  n,  the  electric 
obtained  by  performing  the  operations 


and  magnetic  fields  can  be 


i 


It  follows  from  these  equations  that 


E  =  Idtsinecose  (Y2r2+3yr+3)  e^'r 

^  47r(o+1a)c)  r^ 

V  0  * 


(11.12) 


(11.13) 


(11.14)  _  | 

1 
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E  =  - — -  cos^6(Y^r^+3Yr+3) -(y^r^+yr+1 )  e*^T, 

z  3 

4it(  a+1  u  t )  r  L 

and 

H  =  0  , 

P 

H  =  Milne  (Yr+1)e-vr,  (1 1. 15) 

*  47m 

h2  =  0  • 

where  sin9  »  p/r  and  cosa  =  ( z-£.)/r . 

These  latter  expressions  completely  describe  the  electromagnetic 
fields  produced  in  the  medium  by  a  current  element  source. 
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III.  LINEAR  CURRENT  SOURCE  OF  FINITE  LENGTH 


Suppose  now  that  we  have  an  insulated  thin  straight  wire  extending 
from  to  on  the  2-axis  of  our  coordinate  system  (Figure  III.l)  and 
carrying  a  current  I  exp(iujt),  which  we  assume  is  constant  throughout 
the  length  of  the  wire.  The  surrounding  conducting  medium  is  again 
assumed  to  be  infinite  and  homogeneous  with  electric  and  magnetic  pro¬ 
perties  o,  e,  and  u,  as  described  in  the  previous  section.  We  know  that 
for  a  highly  conducting  medium  such  as  sea  water  (o  =  4  S/m)  it  is  rea¬ 
sonable  to  neglect  the  displacement  currents  for  all  frequencies  up  to 
the  microwave  frequency  range.  We  can  therefore  replace  (o+iuie)  by  a 
in  Equations  (II. 7)  and  (11.11),  obtaining 


Y  - 


(iu >uo  ) 


4 


(III.!) 


"z(c,z) 


jdi 

4^o  r 


(III. 2) 


where,  once  again, 


III.l  Derivation  of  the  Electric  Field  Components 

In  this  subsection  we  follow  the  approach  used  by  Wait  [1952]  to 
derive  the  components  of  the  electric  field  produced  by  the  finite 
linear  current  source.  We  start  with  the  expression  for  the  Hertz  vec¬ 
tor  produced  by  a  current  element  Idf.  (Equation  III. 2).  Integrating 
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from  ^  to  we  obtain  the  Hertz  vector  for  a  finite  linear  current 
source 


n 

z 


(0,2)  » 


4itct 


I (z)e~yr 
r 


d  i  . 


(III. 3) 


Since  we  assumed  that  the  current  is  constant  throughout  the  wire, 
the  Hertz  vector  can  be  rewritten  as 


ilz(p,z) 


(III. 4) 


Substituting  this  in  Equation  (11.12),  it  follows  that  the  p  and  z  com¬ 
ponents  of  the  electric  field  are  given  by 


r 


mri 

3p 


(in. 5) 


E 


z 


r‘2 

I  /  P(r)dt  -  I 
'  *1 


3z 


2 

1  = 


Here  the  function  P(r)  and  Q(r)  are  defined  [Sunde,  1949;  Wait,  1952] 
as  follows: 


p<n  -  ^ 


-yr 


Q(r)  = 


e’2L 

TV  a  r 


(III. 6) 
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Figure  1 1 1 . 1 .  Coordinate  system  and  geometry  used  in  the  derivation  of 
the  electric  anu  magnetic  fields  produced  in  a  conducting  medium  of 
infinite  extent  by  a  linear  current  source  of  finite  extent. 
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Numerical  values  for  these  Integrals  are  tabulated  by  Staff  [1949]. 


The  integral  in  Equation  (III. 5)  can  now  be  written 


slnh'1  j  -  Ec[a,x] 


-  1EsU.x] 


(HI. 9) 


Combining  Equations  (III. 5,  III. 6  and  III. 9),  the  following  final  ex¬ 
pressions  for  the  electric  field  components  are  obtained 


(>'r2+1 ) 


tvr1  +1 ) 


and 


(III. 10) 


r  _  j  vupl 
LZ  4  TT 


!X  X  ) 

sinh'1  ~  -  sinh"1  ~  -  E^a.x.,]  +  Ec[a,x1]| 
-1  |Es[a,x2]  -  Es[a,x1]| 


+  Ti_ 

4ito 


-yr- 


■yr 


S—f-  (yr2+l)  U-*2)  -  — r  (yr^l)  (z-^)  ,  (III. 11) 

r2  rl 


12  2  \  ^ 

where  r^(i=l,2)  is  given  by  r^  =  Ip  +  (z-S,^)  I  . 


1 1 1 . 2  Derivation  of  the  Magnetic  Field  Component 

To  derive  the  magnetic  field  components  produced  by  a  linear 
current  source  of  finite  length>  we  substitute  the  expression  for  the 
Hertz  vector  derived  in  the  previous  section  (Equation  III, 4)  into 


Equation  (11.13).  It  follows  that  the  magnetic  field  has  only  a  single 
component  H  given  by 


y2 

i  14J  4tto 


*2 


L-±  [ 

TTO  3p  J 


d  l  , 


1=1, 


which  can  also  be  written  as 


(III. 12) 


H»  ■  ■  *  L,  [£  ^1 dt  • 


(III  .13) 


Taking  the  derivative  in  the  integral,  we  obtain 


:  4"  J 

r 


H..£f 


(III. 14) 


Making  the  same  substitutions  as  in  the  previous  section,  l.e.,  by 
writing 

6  *  ( 1.11,10/2)**,  y  =  ( 1  +  i )  8.  , 

x  =  6(2-11),  a  =  6p  , 

,2.2' 

u  =  a  +  x 


Equation  (III. 14)  can  be  written  as 


p6  I 

4r 


/. 


? 


-(1  +1  )u  -(l+1)u 


dx  (III. 15) 
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where,  once  again,  x^  a  SU-fc-j)  and  Xj  =  6( 2 - ) .  This  Integral  can  be 
decomposed  Into  the  following  form: 


At  tins  stage,  we  will  define  four  different  Integrals  as  follows: 


x 


(111.16) 
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dx 


e 


u 

cos  u 

T 

u 


g-~U4n-U  dx  . 
u 


Using  these  integrals,  the  magnetic  field  expression  can  be  written  as 


V 


p62I 

4tT 


|Mc[a,  x2]  -  «c[a,  Xl]  +  Ms[a,  x?]  -  Ms[a,  x-j] 
+  Nc[a,  x2]  -  Nc[a,  Xl]j 
♦  1  J«cta,  x2]  -  Mc[a,  x,]  -  M$[a,  x2]  +  M$[a,  x^ 
-  NsCa,  x2]  +  N$[a,  x,]J 


(III. 17) 


with  the  corresponding  component  of  magnetic  induction  B  being  given 

v 

by  the  usual  expression 


V'o"*'  (IUJ8) 

1 1 1. 3  Parametric  Approach 

As  shown  by  Fraser-Smith  and  Bubenlk  [1979],  we  can  make  use  of  a 
parametric  approach  by  normalizing  all  the  distances  In  our  expressions 
by  the  skin  depth  of  the  surrounding  conducting  medium,  6,  where 


20 
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Note  that  the  p^  introduced  here  is  the  same  as  the  quantity  a  used 
previously;  we  introduce  solely  for  notational  consistency. 

These  parametric  equations  (i.e.,  Equat'ions  III. 19,  111.20,  and 
III. 21)  may  be  used  to  convert  the  field  values  calculated  for  one  par¬ 
ticular  frequency  to  another  set  of  field  values  at  some  other  frequency 
without  much  additional  computation.  For  example,  assume  that  the  values 
of  the  electric  and  magnetic  fields  produced  by  the  linear  current  source 


22 


at  the  observation  point  P(p,  <p,  z)  In  Figure  III.l  are  known  for  a 
certain  frequency  f -j .  Let  these  field  values  be  Ep  =  EP1  *  E*  *  Ezl 

and  B.  =  B  -  respectively.  If  we  now  change  the  frequency  by  a  factor 

■P  *P 1 

2  2 

1/A  to  f2  *  f -j /A  the  new  skin  depth  will  be  =  A5-|  .  If  we  multiply 
all  the  distance  coordinates  in  Figure  III.1  by  A,  l.e.,  change  l ^ ,  l*, 
p,  and  z  to  AfE^  ,  ASlj,  Ap,  and  Az,  then  all  the  normalized  distances  in 
the  parametric  equations  stay  the  same  and  the  electric  and  magnetic 
field  components  at  the  new  observation  point  P  (Ap,  <p,  Az)  produced  by 

the  linear  current  source  extending  from  l  -  A«,^  to  l  -  ancl  operating 

2  2  2 
at  the  new  frequency  f^  =  fj/A  become  E  £  =  E^/A  >  E ^  ~  E2i  /A  and 

=  B^/ A.  These  new  field  values  for  frequency  can  obviously  be 

calculated  more  easily  from  the  field  values  for  frequency  f^  than  they 

can  from  the  original  integral  expressions. 

1 1 1 . 4  Computations  of  the  Fields 

The  expressions  derived  in  Sections  III.l  and  III. 2  for  the  field 
components  Ez  and  contain  integrals  which  in  general  can  only  be 
evaluated  numerically.  As  mentioned  in  Section  III.l,  there  are  a  limit¬ 
ed  number  of  tabulated  values  available  [Staff,  1949]  for  the  two  Inte¬ 
grals  E  [a,x]  and  E  [a,x]  in  the  expression  for  E  ,  but  there  are  no 

v  j 

tabulated  values  for  the  four  integrals  Mc[a,x],  Ms[a,x],  Nc[a,x],  and 
Ns[a,x]  in  the  expression  for  B^.  We  therefore  prepared  two  computer 
programs  for  evaluating  (1)  the  two  components  of  the  electric  field, 

Ep  and  Ez>  and  (2)  the  component  of  the  magnetic  field,  B^.  In  these 
programs  the  six  constitutive  integrals  in  the  field  expressions  were 
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evaluated  numerically  by  means  of  Weddle's  rule  [Scarborough,  1966,  p. 
138].  Both  programs  were  interactive  in  the  sense  that  the  value  of  a 
field  component  could  be  evaluated  by  entering  the  location  and  the 
length  of  the  linear  current  source,  the  amplitude  and  the  frequency  of 
the  current  flowing  in  the  source,  and  the  coordinates  of  the  observa¬ 
tion  point  in  the  surrounding  conducting  medium. 

The  values  for  the  numerical  integrals  were  checked  in  several  dif¬ 
ferent  ways.  One  obvious  way  was  by  comparing  our  computed  results  for 
£c[a,x]  and  Es[a,x]  with  the  available  tabulated  values  [Staff,  1949]. 
Another  way  was  to  see  if  the  values  of  the  integrals  converged  to  the 
values  tabulated  for  Kelvin  functions  and  their  derivatives  [Young  and 
Kirk,  1964;  Lowell,  1959;  see  Appendix  B]  for  large  values  of  the  upper 
limit  of  these  integrals.  We  also  used  Simpson's  rule  [Scarborough, 
1966,  p.  137]  and  compared  the  values  obtained  by  two  different  methods 
of  numerical  integration;  namely,  Weddle's  rule  and  Simpson's  rule.  In 
all  cases  our  computed  values  agreed  well  with  the  comparison  values. 

Shown  in  Figures  III. 2  to  III. 7  are  computed  values  of  the  electric 
and  magnetic  fields  produced  at  frequencies  of  100  Hz  (Figures  III. 2, 
III. 4  -  III. 7)  and  1  Hz  (Figures  III. 3,  III. 8)  by  a  linear  current 
source  of  100  m  length  carrying  an  alternating  current  of  1000  A  ampli¬ 
tude.  The  fields  are  shown  in  two  dimensions;  in  all  cases  the  fields 
are  cyl indrical ly  symmetric  about  the  axis  defined  by  the  source.  Fig¬ 
ures  III. 2  and  >11.3  show  the  spatial  variations  of  the  total  electric 
field  vector  at  the  instant  when  there  is  maximum  current  (1000A)  in 
the  source  flowing  from  B  to  A.  The  shorter  wavelength  at  100  Hz  can 
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be  clearly  seen  when  comparison  is  made  of  the  two  figures.  Figure 
1 1 1. 4  shows  equal  amplitude  contours  of  the  magnetic  field  component 

r 

for  a  frequency  of  100  Hz  at  the  time  instant  of  maximum  current  flow 

from  B  to  A.  Figures  III. 5  and  III. 6  show  contours  of  the  maximum 

values  reached  by  the  t*o  components  of  the  electric  field  E_  and  E 

z  P 

during  one  cycle  of  the  source  current.  The  dashed  contours  indicated 
zero  amplitude;  thus  in  Figure  II 1. 6  the  E  component  is  zero  along  the 

p 

axis  of  the  source  and  along  the  perpendicular  line  through  the  center 
of  the  source.  Comparing  the  two  electric  field  figures,  it  is  easy 

to  see  how  the  two  separate  components  E  and  E^  contribute  to  the  total 

P 

electric  field. 

In  Figure  III. 7  we  have  drawn  contours  indicating  the  maximum  values 
reached  by  the  100  Hz  total  magnetic  field  variation  at  each  point  during 

4 

one  cycle.  Thus  at  points  Inside  the  contour  labelled  10  pT,  our  data 
show  that  the  total  magnetic  field  (or,  equivalently,  the  magnetic  field 
component  B  )  will  reach  a  values  of  10^  pT  or  more  during  one  cycle  of 
the  source  current. 

In  Figure  III. 8  we  have  relabelled  tne  data  in  Figure  III. 7  accord¬ 
ing  to  the  parametric  approach  described  in  Section  III. 3;  the  frequency 
has  been  decreased  by  a  factor  of  100,  the  distances  have  been  Increased 
by  a  factor  of  10,  and  the  values  of  on  the  contours  have  been 
decreased  by  a  factor  of  10.  These  changes  leave  the  basic  figure  un¬ 
altered,  but  the  various  contours  now  show  maximum  total  magnetic  field 
values  produced  by  a  source  of  1000  m  length  carrying  a  1000  A  alterna¬ 
ting  current  with  a  frequency  of  1  Hz. 
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Figure  III. 2.  Variation  of  the  total  electric  field  vector  produced  in 
sea  water  by  a  linear  current  source  of  finite  length  (100  m)  carrying 
an  alternating  current  of  amplitude  1000  A  and  frequency  100  Hz  at  the 
time  instant  when  there  is  maximum  current  flowing  in  the  direction  from 
B  to  A.  The  magnitude  of  the  electric  field  vector  is  10n  uv/m,  where  n 
can  be  read  from  the  amplitude  scale  shown  above.  The  pattern  of  elec¬ 
tric  field  vectors  is  cyl indrical ly  symmetric  around  the  axis  of  the 
source . 
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-iqure  III. 3.  Variation  of  the  total  electric  field  vector  produced  in 
sea  water  by  a  linear  current  source  of  finite  length  (100  m)  carrying 
an  alternating  current  of  amplitude  1000  A  and  frequency  1  Hz  at  the 
time  instant  when  there  is  maximum  current  flowing  in  the  direction 
from  B  to  A.  The  magnitude  of  the  electric  field  vector  is  1 0"  uv/m, 
where  n  can  be  read  from  tne  amplitude  scale  shown  above.  The  pattern 
of  electric  field  vectors  is  cyl indrically  symmetric  around  the  axis 
of  the  source. 
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Figure  II I. 4.  Variation  of  the  total  magnetic  field  (which  has  only  the 
one  component  B^)  produced  in  sea  water  by  a  linear  current  source  ot 
finite  length  (100  m)  carrying  an  alternating  current  of  amplitude  1000  A 
and  frequency  100  Hz  at  the  time  instant  when  there  is  maximum  current 
flowing  in  the  direction  from  B  to  A.  The  magnetic  field  magnitudes  are 
given  by  B^  10n  pT,  where  n  is  given  for  each  contour  in  the  figure 
(Ba  =  0  on  the  dashed  contours.)  As  indicated,  the  direction  of  is 
either  into  the  paper  or  out  of  the  paper  depending  on  the  region 
between  the  dashed  lines.  The  field  pattern  is  cyl indrical ly  synmetric 
around  the  axis  of  the  source. 
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Figure  1 1 1. 5.  Variation  of  the  amplitude  of  the  parallel  component  of 
the  electric  field,  Ez>  produced  in  sea  water  by  a  linear  current 
source  of  finite  length  (100  m)  carrying  an  alternating-  current  of 
amplitude  1000  A  and  frequency  100  Hz.  The  amplitude  of  E?  is 
constant  along  each  contour.  The  field  pattern  is  cylindncally  sym 
metric  around  the  axis  of  the  source. 
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Figure  1 1 1. 6.  Variation  of  the  amplitude  of  the  perpendicular  component 
of  the  electric  field,  Er ,  produced  in  sea  water  by  a  linear  current 
source  of  finite  length  '(100  m)  carrying  an  alternating  current  of 
amplitude  1000  A  and  frequency  100  Hz.  The  amplitude  of  E(,  is  constant 
along  each  contour  and  is  zero  along  the  dashed  lines.  The  field  pat¬ 
tern  is  cyl indrical ly  symmetric  around  the  axis  of  the  source. 
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Figure  III. 7.  Variation  of  the  amplitude  of  the  total  magnetic  field 
(which  has  only  the  one  component  B<j,)  produced  in  sea  water  by  a 
linear  current  source  of  finite  length  (100  m)  carrying  an  alternat¬ 
ing  current  of  amplitude  1000  A  and  frequency  100  Hz.  The  field 
pattern  is  cylindrically  symmetric  around  the  axis  of  the  source. 
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Figure  III. 8.  Variation  cf  the  amplitude  of  the  total  magnetic  field 
produced  in  sea  water  by  a  linear  current  source  of  finite  length 
(1000  m)  carrying  an  alternating  current  of  amplitude  1000  A  and 
frequency  1  Hz.  The  field  pattern  is  cyl indrical  ly  symmetric  around 
the  axis  of  the  source  and  is  identical  to  the  pattern  shown  in 
Figure  1 1 1. 7  for  a  frequency  of  100  Hz;  however,  the  scales  have  been 
changed  in  accordance  with  the  rules  for  the  parametric  approach  des¬ 
cribed  in  the  text. 


Finally,  In  Figures  III. 9  and  III. 10  we  have  plotted  data  that  com 
pare  the  electric  and  magnetic  fields  produced  by  the  finite  source  with 
those  produced  by  an  Infinite  source  carrying  the  same  current.  Comparl 
son  of  the  electric  fields  is  made  In  Figure  III. 9,  where  we  show  the 
variation  with  distance  (along  the  perpendicular  axis  through  the  cen¬ 
ter  of  the  finite  source)  of  the  ratio  of  the  amplitudes  of  the  z-com- 
ponents  of  the  electric  field  produced  by  the  finite  (length  100  m)  and 
infinite  sources.  Because  of  the  choice  of  the  central  axis  to  show 
the  distance  variation,  the  Ez  component  for  the  finite  source  is  also 
the  total  electric  field;  we  already  know,  of  course,  that  the  Ez  com¬ 
ponent  for  the  Infinite  source  is  everywhere  the  same  as  the  total  field 
Thus  the  ratio  of  the  Ez  components  shown  in  the  figure  is  also  the 
ratio  of  the  total  electric  fields.  Figure  III. 10  shows  the  variation 
of  the  ratio  of  the  amplitudes  of  the  total  magnetic  field  (i.e.,  Its 
single  component  B^)  for  the  same  two  sources. 

An  interesting  feature  of  the  electric  field  data  in  Figure  III. 9 
is  the  relative  enhancement  of  the  electric  fields  near  the  source  of 
finite  length.  The  increase  only  becomes  apparent,  for  frequencies  less 
than  10  Hz,  but  it  Increases  rapidly  with  decreasing  frequency  and  can 

3 

be  of  the  order  10  at  f  =  0.001  Hz.  The  ratio  of  the  magnetic  fields, 
on  the  other  hand,  never  exceeds  unity  near  the  source  and  it  decreases 
with  increasing  distance  from  the  source. 


Figure  III. 9.  Variation  with  perpendicular  distance  of  the  ratio  of  the 
amplitude  of  the  2-component  of  the  electric  field  produced  by  a 
linear  current  source  of  finite  length  (100  m)  to  the  amplitude  of  the 
z-component  of  the  electric  field  produced  by  a  linear  current  source 
of  infinite  length,  both  carrying  alternating  currents  of  the  same 
amplitude  and  frequency.  The  particular  perpendicular  axis  along 
which  the  variation  is  shown  passes  through  the  center  of  the  finite 
source.  The  surrounding  conducting  medium  (sea  water,  o  =  4  S/m)  Is 
infinite  in  extent.  Each  curve  refers  to  a  different  frequency  In  the 
range  0.001-100  Hz  as  shown  in  the  figure. 
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Figure  III. 10.  Variation  with  perpendicular  distance  of  the  ratio  of 
the  amplitude  of  the  magnetic  field  produced  by  a  linear  current 
source  of  finite  length  (100  m)  to  the  amplitude  of  the  magnetic  field 
produced  by  a  linear  current  source  of  infinite  length,  both  carrying 
alternating  currents  of  the  same  amplitude  and  frequency.  The  parti¬ 
cular  perpendicular  axis  along  which  the  variation  Is  shown  passes 
through  the  center  of  the  finite  source.  The  surrounding  conducting 
medium  (sea  water,  o  =  4  S/m)  is  Infinite  in  extent.  Each  curve 
refers  to  a  different  frequency  in  the  range  0  -  100  Hz  as  shown  in 
the  figure. 
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IV.  LINEAR  CURRENT  SOURCE  OF  INFINITE  LENGTH 


IV. 1  Derivation  of  the  Field  Expressions 

In  the  preceding  section  we  derived  expressions  for  the  electric 
and  magnetic  fields  produced  In  a  conducting  medium  of  Infinite  extent 
by  a  linear  current  source  of  finite  length.  In  the  coordinate  system 
of  Figure  III.l  the  wire  extended  from  2^  to  2^  In  the  z  direction. 
Suppose  now  that  we  let  2^  go  to  and  2j  go  to  +«.  The  electric  field 
components  given  by  Equations  I II. 5  then  become 


Eo  -  0  • 

00 

Ez  -  -I  /  P(r)d2  . 

— ao 


(IV. 1) 


Thus  the  electric  field  of  a  linear  current  source  of  infinite  length 
aligned  in  the  z  direction  has  only  a  z-component.  Making  the  same 
substitution  as  In  Section  III,  the  expression  for  Ez  can  also  be 
written  as 


E,  ■  Wf. 


■(l+1)u 


dx 


2  2  h 

where  u  =  (a  +  x  )  .  Knowing  that  the  Integral  In  this  equation  Is  an 
even  function  of  x,  we  can  also  write 


E 


z 


-1u)U  I  f 
2tt  J  o 


e-(Hi)u 
-  — 


dx  • 


(IV. 2) 


As  described  in  Appendix  8,  the  Integral  can  be  replaced  by  the  modified 


' 
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(IV. 3) 


Bessel  function  of  the  second  kind  of  order  zero,  KQ(yp),  giving 
Ez  =  ^  ko(yp)’ 

which  Is  the  same  expression  as  that  derived  by  Von  Aulock  [1948]  and  by 
Walt  [1959]. 

To  derive  the  magnetic  field,  we  make  use  of  the  above  results  for 
the  electric  field  and  Maxwell's  Equation  (1 1. 3).  After  appropriate 
algebraic  manipulation  we  obtain  the  following  expression  for  the  single 
magnetic  field  intensity  component 

H  ..± 

Differentiation  of  Equation  (iV.3)  then  gives 

V&  »»-4) 

where  K^yp)  is  the  modified  Bessel  function  of  the  second  kind  of  order 
one.  This  latter  expression  for  the  magnetic  field  intensity  was  also 
derived  by  Von  Aulock  [1948]  and  by  Wait  [1959]. 

Summarizing  the  above  results,  we  see  that  a  wire  of  Infinite  length 
carrying  a  uniform  current  I  exp(lwt)  along  the  r  axis  in  an  infinite, 
homogeneous,  conducting  medium  produces  an  electric  field  with  a  single 
component  in  the  z  direction  (i.e.,  directed  parallel  to  the  wire)  given 
by  Equation  (IV.3)  and  a  magnetic  field  with  a  single  compooent  in  the 
<p  direction  (i.e.,  azimuthal  component)  given  by  Equation  (IV. 4). 

IV. 2  Effects  of  Insulation 

In  the  preceding  derivations  of  the  field  expressions  for  linear 
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current  sources  of  finite  and  Infinite  length  we  assume  that  the  current 
is  uniform  throughout  the  length  of  the  conducting  wire  comprising  the 
actual  source  and  that  It  flows  into  the  conducting  medium  only  from  the 
ends  of  the  wire.  For  this  condition  of  uniform  current  to  apply,  it  is 
of  course  necessary  to  have  insulation  on  the  wire.  Without  Insulation 
the  current  in  the  wire  would  flow  into  the  surrounding  conducting  me¬ 
dium  from  all  points  of  the  wire's  surface,  and  uniform  current  flow 
along  the  length  of  the  wire  could  not  be  achieved. 

Extending  results  obtained  by  Sunde  [1949],  Wait  [1952]  states  that 
the  propagation  constant  for  low  frequency  currents  In  an  insulated  wire 
in  a  conducting  medium  is  determined  largely  by  the  electrical  character¬ 
istics  of  the  Insulation,  and  he  gives  the  following  approximate  expres¬ 
sion  for  the  current  I ( x )  at  a  point  In  the  wire  located  a  distance  x  away 
from  the  generator  terminals: 

I(x)  ~  IQ  exp  (-rx), 

where  IQ  Is  the  current  assumed  to  be  entering  the  wire.  The  propagation 
constant  Tin  this  expression  is  given  approximately  by 

r  2  lU^u2)5*, 

where  and  are  the  permittivity  and  permeability  of  the  Insulation 

material.  It  follows  from  these  expressions  that  the  current  along  an 

Insulated  wire  of  length  L  will  be  essentially  uniform  provided  JTL(  «  i. 

Suppose  we  take  polyethylene  to  be  a  typical  insulating  material.  For 

this  material  we  can  write  e.  =  e  e  .  with  e  =  2.25  (the  dielectric  con- 

i  r  o  r 

stant  of  polyethylene),  and  p  *  uQ.  Substituting  the  known  values  of 
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eQ  and  pQ,  we  obtain 

| r |  -  (3.14  x  10'8f)  m'1  , 

-8  -1 

For  f  -  1  Hz,  |r|  *  3  x  10  m  and  the  length  of  wire  can  be  as 
large  as  1000  km  and  yet  the  condition  J TL |  «  1  is  still  satisfied. 
Thus  for  low  frequencies  and  typical  insulating  materials,  it  is  reason¬ 
able  to  assume  that  an  alternating  current  is  uniform  throughout  the 
length  of  a  wire  for  wire  lengths  varying  from  several  hundreds  to 
several  thousands  of  kilometers. 

It  also  follows  that  the  displacement  currents  flowing  to  the 
surrounding  medium  through  the  insulation  are  negligible  and  the  current 
mostly  flows  to  tht  medium  from  the  ends  of  the  wire. 

It  is  also  shown  analytically  by  Wait  [1952]  and  stated  by  Kraich- 
man  [1976]  that  the  electric  and  magnetic  field  expressions  for  linear 
current  sources  of  both  finite  and  infinite  length  (with  the  wire 
assumed  to  be  of  negligible  thickness  in  both  cases)  are  unaffected  by 

9 

the  properties  of  the  insulating  material  on  the  wire  as  long  as  the 
ratio  of  the  radius  of  the  insulation  to  the  wavelength  of  the  surround¬ 
ing  conducting  medium  is  much  less  than  one.  This  condition  can  be 
written  in  the  form 

b(fyo)15  «  1 , 

where  b  is  the  radius  of  insulation  and  p  and  a  are  the  permeability  and 
the  conductivity  of  the  surrounding  medium.  For  sea  water  (o  =  4  S/m, 
p  =  pQ  =  4ir  x  lO"7  H/m)  this  condition  becomes 
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2.24  x  10-3  bflj«l, 

which  will  of  course  be  satisfied  for  low  frequencies  and  reasonable 
thickness  of  insulation. 


IV. 3  Computation  of  the  Fields 

Plots  of  the  variations  with  distance  of  the  amplitudes  of  Ez  and 
produced  in  sea  water  by  linear  current  sources  of  Infinite  length 
are  given  by  Von  Aulock  [1948]  and  Kraichman  [1976]  for  a  few  selected 
frequencies  and  for  limited  ranges  of  magnitude  and  distance.  These 
data  play  such  an  Important  role  in  the  study  of  the  electric  and  mag¬ 
netic  fields  produced  In  sea  water  by  submerged  linear  current  sources 
that  we  have  extended  the  Von  Aulock/Kraichman  curves  over  a  number  of 
decades  In  frequency,  amplitude  and  distance. 

In  our  derivation  of  the  numerical  data  for  these  curves,  we  made 
use  of  the  fact  that  Equations  (IV. 3)  and  (IV. 4)  can  be  rewritten  In 
terms  of  the  Kelvin  functions  and  their  derivations  to  give  the  follow¬ 
ing  two  new  equations 


r  _  1 U)  U  I 

Ez  '  TT~ 


kerQ  (v^a)  +  1  keiQ  [\fz  a) 


(IV. 5) 


and 


H  -  -  -M. 

♦  wl 


ker^  (\/z  a)  +1  keV  ( \/2~  a )  J  ,  (IV. 


6) 
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W 

where  a  =  up  and  p  =  (wuc/2)  ,  as  before.  It  follows  that  Ez  and 
can  be  evaluated  either  by  using  tables  of  the  Kelvin  functions  and 
their  derivatives  [Lowell,  1959;  Young  and  Kirk,  1964]  in  conjunction 
with  Equations  (IV. 5)  and  (IV. 6),  or  by  using  numerical  integration 
techniques.  We  used  both  methods  to  obtain  our  numerical  values  for 
the  amplitudes  of  E  and  H  . 

Z  9 

The  curves  we  obtained  for  the  variation  of  the  amplitude  of  the 
electrical  field  component  Ez  with  frequency  and  perpendicular  distance 
from  the  wire  are  shown  in  Figure  IV. 1,  and  the  corresponding  curves 
for  the  amplitude  of  the  magnetic  field  component  are  shown  in  Fig¬ 
ure  IV. 2.  In  these  curves  the  lower  limits  of  the  amplitudes,  10"^  mV/m 
for  the  electrical  field  and  10'^  pT  for  the  magnetic  field,  are  roughly 
the  smallest  values  that  can  be  measured  at  depth  in  the  sea  with  pre¬ 
sent  technology. 

There  is  an  interesting  feature  of  the  electric  field  variation 
that  particularly  distinguishes  it  from  the  magnetic  field  variation: 
the  electric  field  first  increases,  reaches  a  maximum  value,  and  then 
decreases  with  increasing  frequency,  whereas  the  magnetic  field  decreases 
monotonical ly  with  increasing  frequency.  This  feature  can  be  seen 
more  clearly  when  the  variation  of  the  amplitude  of  the  electric  field 
Is  plotted  against  frequency,  as  shown  in  Figure  IV. 3,  and  it  implies 
that,  for  a  given  perpendicular  distance  between  the  source  and  the 
field  point,  there  is  an  optimum  frequency  at  which  the  amplitude  of 
the  electric  field  is  a  maximum  at  the  field  point. 


42 


0.001  0.01  0.1  1  10 

DISTANCE  (km) 


Figure  IV.  1.  Variation  with  perpendicular  distance  of  the  amplitude  of 
the  electric  field  produced  in  sea  water  by  a  linear  current  source 
of  infinite  length  carrying  an  alternating  current  of  amplitude  1000A 
The  electric  field  amplitudes  for  an  arbitrary  current  I  can  be 
derived  from  the  values  given  above  by  multiplying  by  I/1Q00. 
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Figure  IV. 2.  Variation  with  perpendicular  distance  of  the  amplitude  of 
the  magnetic  field  produced  in  sea  water  by  a  linear  current  source 
of  infinite  length  carrying  an  alternating  current  of  amplitude  1000  A. 
The  magnetic  field  amplitudes  for  an  arbitrary  current  I  can  be 
derived  from  the  values  given  above  by  multiplying  by  1/1000. 
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Figure  IV. 3.  Variation  with  frequency  of  the  amplitude  of  the  electric 
field  produced  in  sea  water  by  a  linear  current  source  of  infinite 
length  carrying  an  alternating  current  of  amplitude  1000  A  at  a  fixed 
perpendicular  distance  from  the  source.  The  electric  field  amplitudes 
for  an  arbitrary  current  can  be  derived  from  the  values  given  above  by 
multiplying  by  1/1000. 


IV. 4  Parametric  Representation  of  the  Fields 


It  Is  possible  to  write  the  amplitudes  of  the  electric  and  magne¬ 
tic  fields  given  by  Equations  (IV. 5)  and  (IV. 6),  in  the  parametric  form, 
using  the  approach  discussed  in  Section  III. 3.  Normalizing  the  perpen¬ 
dicular  distance  of  the  field  points  from  the  sources  by  the  skin  depth, 
and  generally  following  the  rules  in  Section  III. 3,  the  electric  and 
magnetic  ield  expressions  can  be  written  in  the  following  parametric 
form  : 

{  it  a  5  2/I )  J  Ez|  =  (ker2  (  Jl  a  )  +  kei2  (JT  tft  %  (IV. 7) 

and 

(v6-5/u0I )  I  B^'  =  (ker^2  (Jit)  +  kei’2  (Jit)}  .  (IV. 8) 

The  right  hand  sides  of  these  two  parametric  expressions  can  be 
evaluated  without  difficulty  and  the  data  that  are  obtained  are  plotted 
in  Figure  IV. 4.  The  curve  with  crossed  circles  gives  the  variation  of 
the  amplitude  of  the  electric  field  in  parametric  form,  as  given  by 
Equation  (IV. 7)  and  the  other  curve  (dotted  circles)  gives  the  varia¬ 
tion  of  the  amplitude  of  the  magnetic  field  In  parametric  form,  as 
given  b>  Equation  (IV. 8).  The  variable  along  the  horizontal  axis  In 
these  figures,  \Jl  a ,  is  the  argument  of  the  Kelvin  functions;  It  Is 

essentially  the  perpendicular  distance  from  the  source  in  terms  of  skin 

depths.  The  two  curves  if.  Figure  IV. 4  are  useful  for  rapid  computation 
of  the  amplitudes  of  either  the  electric  or  magnetic  fields  produced  in 
any  arbitrary  conducting  medium  of  infinite  extent  by  a  linear  current 
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Fiqure  IV. 4.  Variation  with  parametric  distance  Jl  a  (where  a  -  P<5  - 
p/6)  of  the  parametric  expressions  for  the  amplitude  of  the  electric 
field  (crossed  circles)  and  the  magnetic  field  (dotted  circles)  Pro¬ 
duced  in  a  conducting  medium  of  infinite  extent  by  a  1  inear  current 
source  of  infinite  length  carrying  an  alternating  current  of  amplitude 

t. 
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source  of  infinite  length  over  a  very  wide  range  of  source  currents 

and  frequencies.  For  example,  if  the  observation  point  Is  100  m  away 

(perpendicular  distance)  from  the  source  and  the  source  frequency  is 

1  Hz,  we  have^/Ta  =  0.56  and  using  Figure  IV. 4  we  obtain  approximately 

1.0  and  1.6  for  the  parametric  electric  and  magnetic  field  quantities 

given  on  the  vertical  axis  of  the  figure.  Knowing  that  6  =  251.6  m  for 

-7  3 

f  =  1  Hz,  a  =  4  S/m,  and  yQ  =  4tt  x  10  H/m,  and  using  I  =  10  A,  we 
obtain  |E2|  =  1.26  mV/m  and  |B^|  =  1 .80  x  10^  pT.  These  particular 
values  can  be  checked  against  the  values  of  | |  and  |B^|  plotted  in 
Figures  IV. 1  and  IV. 2. 

We  have  also  plotted  in  Figure  IV. 5  the  amplitude  of  the  electric 
field  in  the  alternate  parametric  form 

2 

(2^ao2/I)  |EZI  *  (  y/T  a)  (  ker2  ( /T  a)  +  kei2(  /fa)  )  .  (IV. 9) 

This  is  basically  the  same  equation  as  Equation  (IV. 7),  except  that  in 
Equation  (IV. 7)  the  left  hand  side  is  not  a  function  of  the  perpendicu¬ 
lar  distance  p  and  the  left  hand  side  of  the  above  equation  is  not  a 
function  of  frequency.  In  other  words,  the  single  curve  shown  in  Fig¬ 
ure  IV. 5  is  a  parametric  form  of  the  curves  seen  in  Figure  IV. 3;  it 
contains  all  their  information  on  the  variation  of  the  amplitude  of  the 
electric  field  with  perpendicular  distance  and  source  frequency,  as 
well  as  much  additional  information  of  the  field  amplitudes  at  other 
distances.  (If  required,  the  field  amplitudes  can  also  be  computed  in 
other  media  than  sea  water.)  The  maximum  value  of  the  parametric  electric 
field  expression  occurs  for  JT  a  -  2.17.  This  is  a  useful  result  to 
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Figure  IV. 5.  Variation  with  parametric  distance  a  (where  a  =  o,  = 

2  u 
p/6)  of  the  parametric  expression  2nop  |E  |/I  for  the  amplitude  of 

the  electric  field  produced  In  a  conducting  medium  of  Infinite  extent 

by  a  linear  current  source  of  infinite  length  carrying  an  alternating 

current  of  amplitude  I.  Note  that  the  maximum  parametric  amplitude 

(0.83)  occurs  for  JY  a  =  2.17. 


remember,  because  it  enables  the  optimum  frequency  of  operation  mentioned 
in  the  previous  section  and  the  corresponding  maximum  electric  field 
to  be  computed.  The  procedure  to  obtain  this  information  Is  as  follows. 
First,  knowing  the  location  of  the  observation  point  (i.e.,  the  perpen¬ 
dicular  distance  p),  It  Is  simple  to  obtain  the  frequency  maximizing  the 
amplitude  of  the  electric  field  at  that  specific  point  from  JY  a  * 

J 2  p/6  =  2.17.  This  is  the  frequency  we  will  refer  to  as  the  opti¬ 
mum  frequency.  Next,  the  data  in  Figure  IV. 5  show  that  the  maximum 
value  of  the  parametric  electric  field  expression  2irap  |  Ez|  /I  is  approx¬ 
imately  0.83,  from  which  the  actual  value  of  the  electric  field  ampli¬ 
tude  can  be  computed. 


V.  LINEAR  CURRENT  SOURCE  OF  SEMI -INFINITE  LENGTH 


V.1  Derivation  of  the  Electric  and  Magnetic  Field  Components 

In  the  previous  section  we  obtained  expressions  for  the  electric 
and  magnetic  fields  produced  by  a  linear  current  source  of  infinite 
length  by  letting  2,^  and  tend  to  minus  Infinity  and  Infinity,  respec¬ 
tively,  In  the  expressions  for  the  fields  produced  by  a  linear  current 
source  of  finite  length.  In  this  section  we  concentrate  our  attention 
on  the  electric  and  magnetic  fields  surrounding  the  end  of  a  linear 
current  source  Immersed  in  a  conducting  medium  of  Infinite  extent,  l.e., 
around  an  electrode  of  negligible  size  In  the  medium.  We  do  this  by 
letting  just  one  end  of  the  finite  linear  current  source  tend  to  infi¬ 
nity,  thus  producing  a  linear  current  source  of  semi-infinite  length. 

Starting  with  the  finite  source  electric  and  magnetic  fields  ex¬ 
pressions  (Equations  III. 5  and  III. 12),  letting  2^  go  to  -®,  and  gen¬ 
erally  following  a  similar  approach  to  the  one  used  In  Sections  III  and 
IV,  we  obtain  the  following  three  expressions  for  the  non-zero  fields 
In  the  conducting  medium: 

EP  ♦  >>•  l*’’) 

r2 


1  UJU  I 

Tir 


jslnh"1  ~  -  Ec  [a,  x^]  -  kerQ 


(</*■•)! 


-1  j  Es  u,  X2]  +  kel 0  ( sfl  a)j 
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+  ^  ^~~r  {yr2 + 1}  (z  •  V- 


(V.2) 


H  -  -  fif- 
4>  4tt 


1  jMcCa* 


x2]  +  M$[a ,  x2]  +  Nc[a,  x2] 


+i  |MCU,  x2]  -  Ms[a,  x2]  -  Ns[a,  x2] 


-tt- fker0‘  (V?  a)  +  ikey  (Vfa)  , 

zjl  rr  J 


(V.  3) 


where  all  the  terms  here  are  defined  as  In  Section  III. 

When  the  observation  point  Is  on  the  axis  ot  the  semi -Infinite 
wire  (l.e.,  p  =0  and  2  >H2),  the  above  expressions  simplify  greatly  to 


Ep-°- 

Ez  =  ‘1^  jve  ♦  7  1n2  +  lnx2  -  Ecto,  x2]j  -1  { E$[0,  x2]  -  J-  J 
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In  these  last  equations  Ye  Is  Euler's  constant,  which,  to  four  signifi¬ 
cant  figures,  is  given  by  yg  =  0.5772.  This  constant  Is  usually  denoted 
by  the  symbol  y,  but  we  have  added  the  subscript  e  to  distinguish  it 
from  the  y  used  here  for  the  propagation  constant.  There  is,  of  course, 
no  relation  between  Euler's  constant  yg  and  the  propagation  constant  y. 

V.2  Some  Numerical  Values 

Figure  V.l  shows  contours  of  the  maximum  values  reached  by  the  two 
components  of  the  electric  field,  Cp  (solid  line)  and  E2  (dashed  line) 
around  the  end  A  of  the  semi -inf inite  source  during  one  cycle  of  the 
current  in  the  source  (the  amplitude  of  the  current  is  assumed  to  be 
1000  A).  As  can  be  seen,  there  are  some  regions  in  the  figure  where 
the  p-component  is  larger  than  the  z-component  of  the  electric  field, 
but  for  large  distances  the  z-component  becomes  by  far  the  largest  com¬ 
ponent.  This  difference  at  large  distances  Is  caused  by  the  overall 
inverse  cube  decline  of  the  p-component  with  distance  (see  Equation  V.l), 
in  contrast  to  the  z-component,  which  has  only  a  partial  Inverse  cube 
dependence  on  distance  (Equation  V.2).  Note  that  as  the  point  of  obser¬ 
vation  begins  to  move  towards  the  other  end  of  the  source,  i.e.,  the  end  at 
Infinity,  the  z-component  contours  become  parallel  to  the  source  and  the 
electric  field  values  tend  toward  those  produced  by  a  linear  current 
source  of  infinite  length. 
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Figure  V.l.  Variation  of  the  amplitudes  of  the  p  and  z  components  of 
electric  field,  Eft  and  E,,  produced  in  sea  water  by  a  linear  current 
source  of  semi -infinite  length  (extending  from  to  A)  and  carrying 
an  alternating  current  of  amplitude  1000  A  and  frequency  100  Hz.  The 
solid  contours  in  this  figure  correspond  to  Ep,  the  dashed  contours 
correspond  to  E2,  and  the  values  of  n  shown  for  each  contour  give  the 
numerical  amplitude  of  the  corresponding  electric  field  in  units  of 
10n  mV/m.  Note  that  the  field  patterns  are  cylindrically  symmetric 
about  the  axis  defined  by  the  source. 


VI.  LINEAR  CURRENT  SOURCE  OF  INFINITE  LENGTH  WITH  A  FINITE  GAP 


VI .1  Derivation  of  the  Fields 

In  Section  IV  we  investigated  the  electromagnetic  field  produced 
by  an  alternating  current  flowing  in  an  insulated  continuous  linear 
current  source  of  infinite  length  immersed  in  a  homogeneous  conducting 
medium  of  infinite  extent.  This  system  is  distinguished  by  the  fact 
that  all  the  impressed  current  is  effectively  confined  to  the  source. 

In  other  systems,  for  example  linear  current  sources  of  finite  or  semi¬ 
infinite  length,  the  impressed  current  also  flows  into  the  surrounding 
conducting  medium  from  the  ends  of  the  wire.  Suppose  we  now  introduce 
a  gap  into  our  linear  current  source  of  infinite  length  so  that  the 
current  in  the  source  is  forced  to  flow  through  the  conducting  medium 
over  some  finite  distance.  This  system  is  equivalent  to  two  semi- 
infinite  linear  current  sources  carrying  the  same  current  and  aligned 
along  the  same  axis  but  with  a  finite  gap  between  them  as  shown  in 
Figure  VI. 1(a). 

We  have  already  derived  the  electric  and  magnetic  field  expressions 
for  a  semi-infinite  linear  current  source,  so,  using  superposition,  we 
can  easily  show  that  the  fields  produced  by  the  system  in  Figure  VI. 1(a) 
are  given  by  the  following  equations: 
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Figure  VI. 1.  Sketches  of  the  two  linear  current  sources  of  infinite 
length  with  gaps  that  are  considered  in  this  work.  In  each  case  the 
gap  extends  from  l  =  to  l  *  ^  and  the  amplitude  of  Uc  current  is 
the  same  in  both  the  upper  and  lower  semi-infinite  segments.  However, 
the  directions  of  the  current  flow  are  different  in  the  two  systems. 

In  case  (a)  the  currents  flowing  In  each  segment  are  In  the  same 
direction,  whereas  in  case  (b)  the  currents  are  oppositely  directed. 
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It  Is  Interesting  to  note  that,  once  again  as  a  result  of  the  prin¬ 
ciple  of  superposition,  these  expressions  can  be  derived  by  subtracting 
the  field  expressions  for  a  finite  source  of  the  same  length  as  the  gap 
from  the  field  expressions  for  the  continuous  source  of  infinite  length. 
Furthermore,  if  we  are  on  the  axis  of  the  source  and  within  the  gap. 
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(i.e.,  p  =  o  and  the  field  expressions  become 


A  different  version  of  this  gap  system  can  be  obtained  simply  by 
reversing  the  direction  of  the  current  in  one  of  the  semi-infinite  seg¬ 
ments,  so  that  the  currents  are  oppositely  directed  in  the  two  segments 
(Figure  VI. 1(b)).  For  this  case,  we  can  again  make  use  of  the  semi¬ 
infinite  wire  expressions  to  aerive  the  following  electric  and  magnetic 
field  expressions: 
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As  we  will  see  in  the  following  section,  the  fields  produced  by  the  two 
gap  systems  have  interesting  differences  in  their  properties. 


VI .2  Numerical  Results 

As  can  be  seen  from  the  form  of  the  field  expressions  VI. 1  -  VI. 3 
and  VI. 5  -  VI. 7,  the  spatial  variations  of  the  fields  produced  in  and 
around  the  two  finite  gaps  considered  in  this  work  are  not  simple.  To 
illustrate  these  spatial  variations,  and  to  provide  some  representative 
numerical  data,  we  make  the  comparatively  simple  fields  produced  by  a 
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continuous  infinite  linear  current  source  our  reference  and  we  compare 
the  gap  fields  with  these  reference  fields  by  computing  ratios  of  the 
corresponding  field  quantities. 

The  first  of  these  comparisons  is  made  in  Figure  VI. 2.  where  we 
have  plotted  the  variation  with  perpendicular  distance  of  the  ratio  of 
the  ampl;tudes  of  the  total  electric  field  produced  by  (a)  the  infinite 
linear  current  source  with  a  gap  illustrated  in  Figure  VI. 1(a)  (current 
flow  ng  in  the  same  direction  in  the  two  semi-infinite  segments)  and  by 
(b)  the  continuous  infinite  linear  current  source.  The  current  is  every¬ 
where  assumed  to  be  the  same  in  the  sources,  the  gap  width  is  taken  to 
be  100  m,  and  the  variation  with  distance  is  taken  along  the  perpendicu¬ 
lar  axis  passing  through  the  center  of  the  gap.  Because  ui  Lb"  choice 
of  this  particular  axis  to  show  the  distance  variation,  the  p-component 
of  the  electric  'ield  produced  by  the  gap  system  is  zero  at  all  the  field 
points  that  are  considered  and  the  total  electric  field  and  Ez  are  equi¬ 
valent  for  both  sources. 

The  most  interesting  feature  of  the  data  in  Figure  VI. 2  is  the 
large  ratio  at  small  distances  (i.e,,  p<  1000  m)  for  f  <  1  Hz.  At  these 
frequencies  the  amplitude  of  the  electric  field  in  the  vicinity  of  the 
gap  can  be  substantially  larger  than  the  amplitude  of  the  field  that 
would  be  produced  at  the  same  locations  by  a  continuous  source  carrying 
the  same  current.  The  increase  begins  to  decline  rapidly  once  a  minimum 
perpendicular  distance  is  exceeded.  In  Figure  VI. 2  this  distance  is 
aporoximately  50  m,  which  also  happens  to  be  the  distance  from  the  cen¬ 
ter  of  the  gap  to  the  r.ear  ends  of  the  two  semi-infinite  segments. 
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Figure  VI. 2.  Variation  with  perpendicular  distance  of  the  ratio  of  the 
amplitudes  of  the  z-components  of  the  electric  field  produced  (a)  by 
a  linear  current  source  of  infinite  length  with  a  finite  gap  (100  m), 
carrying  current  flowing  in  the  same  direction  in  each  semi -infini te 
segment,  and  (b)  by  a  continuous  linear  current  source  of  infinite 
length.  Both  sources  carry  alternating  currents  of  the  same  amplitude 
and  frequency  and  the  perpendicular  axis  along  which  the  distance 
variation  is  shown  passes  through  the  center  of  the  gap.  The  surround¬ 
ing  conducting  medium  (sea  water,  a  =  4  S/m)  is  assumed  to  be  infinite 
in  extent. 
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Beyond  this  minimum  distance  the  ratio  of  the  total  field  amplitudes 
rapidly  approaches  unity,  no  matter  what  frequency  is  involved,  and  at 
large  distances  the  fields  of  the  gap  system  approximate  those  produced 

by  the  continuous  linear  current  source  of  infinite  length. 

A  second  comparison  of  the  gap  fields  with  the  reference  fields  pro¬ 
duced  ; y  a  continuous  source  is  made  in  Figure  VI. 3.  In  this  figure  the 
gap  is  the  same  size  (100  m)  as  in  the  previous  example,  but  the  currents 
in  the  two  semi-infinite  segments  are  now  taken  in  opposite  directions. 
The  figure  therefore  shows  the  variation  with  perpendicular  distance  of 
the  ratio  of  the  amplitudes  of  the  total  electric  field  produced  by  (a) 
the  infinite  linear  current  source  with  a  gap  illustrated  in  Figure  IV. 
1(b)  (current  flowing  in  opposite  directions  in  the  two  semi-infinite 
segments)  and  by  (b)  the  continuous  infinite  linear  current  source.  As 
before,  the  amplitude  of  the  current  is  everywhere  considered  to  be  the 
same  in  the  sources  and  the  variation  with  distance  is  taken  along  the 
perpendicular  axis  passing  through  the  center  of  the  gap.  Again  because 
of  the  choice  of  this  particular  axis  to  show  the  distance  variation, 
the  electric  fields  produced  by  the  gap  system  have  a  distinguishing 
characteristic:  the  z-component  of  the  electric  field  is  zero  at  all 
the  field  points  that  are  considered  and  the  total  electric  field  and 
are  equivalent.  For  the  continuous  source,  the  total  electric  field 
and  ^z  continue  to  be  equivalent. 

As  was  observed  in  the  previous  example,  there  is  a  region  near  the 
gap  where  the  ratio  is  much  larger  than  one  for  the  lower  frequencies 
(f  <  1  Hz).  This  region  extends  out  somewhat  further  from  the  gap  than 
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Fig.  VI. 3.  Variation  with  perpendicular  distance  of  the  ratio  of  the 
amplitudes  of  (a)  the  p-component  of  the  electric  field  produced  by  a 
linear  length  with  a  finite  gap  (100  m),  with  current  flowing  in  oppo¬ 
site  direction  in  each  semi-infinite  segment,  and  (b)  the  z-component 
of  the  electric  field  produced  by  a  continuous  linear  current  source 
of  infinite  length.  Both  sources  carry  alternating  current  of  the 
same  amplitude  and  frequency  and  the  perpendicular  axis  along  which 
the  distance  variation  isshown  passes  through  the  center  of  the  gap. 
The  surrounding  conducting  medium  (sea  water,  a  =  4  S/m)  is  assumed 
to  be  infinite  in  extent. 
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was  observed  before,  which  may  be  considered  a  positive  characteristic  of 
this  particular  gap  system’s  electric  field.  On  the  other  hand,  at  large 
distances  the  gap  system's  electric  field  decays  more  rapidly  with  dis¬ 
tance  than  the  electric  field  of  the  reference  source  and  the  ratio 
plotted  in  the  figure  becomes  less  than  one. 

Our  third  and  final  comparison  of  the  electric  fields  produced  by  a 
gap  system  with  the  fields  produced  by  other  systems  is  made  in  Figure 
VI. 4.  In  this  case  we  have  plotted  the  variations  with  frequency  of  the 
total  electric  fields  produced  at  three  different  observation  points,  as 
follows:  (1)  an  observation  point  at  the  center  of  the  gap  system  (gap 
width  -  100  m)  illustrated  in  Figure  VI. 1(a);  (2)  an  observation  point 
50  m  away  from  the  center  of  a  linear  current  source  (length  «  100  m) 
and  located  on  a  perpendicular  line  through  the  center  of  the  source; 
and  (3)  an  observation  point  located  50  m  away  from  a  continuous  linear 
current  source  of  infinite  length  (our  reference  source).  In  all  cases 
the  amplitude  of  the  current  in  the  sources  is  1000  A. 

These  final  data  show  once  again  how  the  electric  field  of  the 
reference  source  can  be  exceeded  substantially  by  the  fields  produced 
either  by  gap  systems  or  sources  of  finite  length  carrying  the  same  total 
current  when  the  frequencies  and  distances  are  small.  However,  as  the 
field  point  distances  Increase  and  become  large  in  terms  of  skin  depths 
the  fields  produced  by  the  continuous  linear  current  source  of  infinite 
extent  are  at  least  as  large  as  any  of  the  other  fields,  and  they  may  be 
substantially  larger  in  some  cases. 
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Figure  VI. 4.  Comparison  of  the  variations  with  frequency  of  the  ampli¬ 
tudes  of  the  z-component  of  the  electric  field  for  three  cases:  (1) 

}E2j  produced  by  a  linear  current  source  of  Infinite  length  with  a 
finite  gap  (100  m)  carrying  alternating  current  of  amplitude  1000  A 
and  with  the  observation  point  fixed  at  the  center  of  the  gap  (dotted 
circles),  (2)  jEz|  produced  by  a  linear  current  source  of  finite 
length  (100  m)  carrying  an  alternating  current  of  amplitude  1000  A, 
with  the  observation  point  fixed  50  m  along  a  perpendicular  line  through 
the  center  of  the  source  (triangles),  (3)  |EZ|  produced  by  a  continuous 
linear  current  source  of  infinite  length  carrying  an  alternating  1000  A, 
with  the  observation  point  fixed  50  m  away  from  the  source  (crossed 
circles).  The  surrounding  conducting  medium  (sea  water,  o  =  4  S/m) 
is  assumed  to  be  infinite  in  extent.  For  each  of  the  three  cases  tne 
amplitude  of  the  total  electric  field  and  |£ZJ  are  equivalent. 
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The  electric  and  magnetic  fields  along  the  straight  line  joining 
the  two  ends  of  the  gap,  i.e. ,  along  the  axis  of  the  gap,  have  a  simple 
vector  form:  the  electric  field  has  only  a  single  component  directed 
along  the  axis  and  the  magnetic  field  is  zero.  Unfortunately,  this 
simplicity  does  not  carry  over  to  the  spatial  variation  of  the  electric 
field  along  the  axis.  Further,  just  off  the  axis  the  magnetic  field 
has  a  non-zero  azimuthal  component,  and  this  component  also  has  a  com¬ 
plicated  spatial  variation  as  the  observation  point  moves  along  a  line 
parallel  to  the  axis.  In  general,  both  the  axial  electric  field  compo¬ 
nent  and  the  off-axis  magnetic  field  component  are  large  near  the  two 
ends  of  the  gap  and  they  decrease  as  the  observation  point  moves  toward 
the  center  of  the  gap.  The  magnitude  of  the  decrease  depends  in  parti¬ 
cular  on  the  width  of  the  gap  and  on  the  frequency:  if  the  gap  width  is 
less  than  a  skin  depth,  the  decrease  can  be  quite  small;  but  If  the 
width  exceeds  a  skin  depth,  the  decrease  can  be  very  large  near  the 
center  of  the  gap.  This  behavior  is  well  illustrated  by  the  curve  in 
Figure  VI. 4  showing  the  variation  with  frequency  of  the  amplitude  of 
the  total  electric  field  (I.e.,  the  axial  or  z-component)  at  the  center 
of  a  100  m  gap  in  sea  water.  It  can  be  seen  that  as  the  frequency  in¬ 
creases  from  0.001  Hz  there  is  little  variation  In  the  amplitude  until 
the  frequency  reaches  about  4  Hz  (skin  depth  =  125.8  m),  but  by  the  time 
the  frequency  has  reached  10  Hz  (skin  depth  =  79.6  m)  the  amplitude  has 
noticeably  begun  to  decline. 
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VII.  LINEAR  CURRENT  SOURCE  OF  INFINITE  LENGTH  WITH  AN  ELEMENTARY  GAP 


In  the  previous  section  we  derived  the  fields  produced  by  two  semi¬ 
infinite  linear  current  sources  aligned  along  the  same  axis  with  a  gap 
between  them.  If  we  now  reduce  the  size  of  the  gap  until  it  becomes  in¬ 
finitesimal,  we  will  obtain  a  linear  current  source  of  infinite  length 
with  an  elementary  gap  in  the  middle. 

Suppose  the  gap  has  an  elementary  length  x  and  that  it  is  centered 
on  the  orgin.  Starting  from  the  Hertz  vector  expression  for  a  current 
element  source  (Equation  III. 2)  and  integrating  it  with  the  right  limits, 
we  obtain  the  total  Hertz  vector  due  to  the  infinite  length  wire  with  an 
elementary  gap: 


which  can  also  be  written  in  the  form 


T<»  T  /  C. 

,1/*  di .  i  c  dL 

4tto  j  r  '4xo  J  ._  r 


(VII. 2) 


The  two  terms  on  the  right  of  this  last  expression  are  easily  recog¬ 
nized:  the  first  term  is  the  Hertz  vector  due  to  a  continuous  infinite 
linear  current  source  and  the  second  term  is  the  Hertz  vector  due  to  a 
current  element  of  elementary  length  x. 

Using  Equations  (II. 5)  and  (II. 6),  the  electric  and  magnetic  fields 
tor  this  system  can  be  derived  with  little  difficulty;  they  are  given  by 
the  following  expressions:  • 
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where  sine  =  p/r,  cose  =  z/r,  and  r  =  (p  +  z  )  .  It  is  obvious,  of 
course,  that  these  field  expressions  are  simply  those  for  a  continuous 
infinite  wire  less  the  corresponding  expressions  for  a  current  element 
of  length  t  centered  on  the  origin. 


VIII.  SUMMARY  AND  CONCLUSIONS 


VIII. 1  Summary 

Starting  with  the  theoretical  work  of  Von  Aulock  [1948,  1953],  Wait 
[1952,  1959,  1960],  and  Kraichman  [1976],  we  have  derived  a  number  of 
new  theoretical  and  numerical  results  for  the  ULF/ELF  electromagnetic 
fields  produced  in  a  conducting  medium  of  infinite  extent  by  a  variety 
of  linear  current  sources.  Many  of  the  numerical  results  apply  specifi¬ 
cally  to  a  sea  water  medium,  since  our  work  is  oriented  toward  an  under¬ 
sea  communication  application,  but  the  basic  theory  and  some  of  the 
numerical  data  for  the  parametric  expressions  are  given  in  a  form  that 
is  independent  of  the  specific  conducting  medium.  The  following  is  a 
list  of  the  new  results  we  have  derived;  in  all  cases  the  sources  and 
observation  points  are  taken  to  be  in  a  conducting  medium  of  infinite 
extent: 

(1)  We  have  greatly  expanded  the  available  numerical  data  on  the 
electrical  and  magnetic  field  amplitudes  produced  by  linear  current 
sources  of  infinite  length  for  frequencies  in  the  range  0.0C1  -  1000  Hz. 

We  have  also  introduced  a  parametric,  representation  for  these  fields  and, 
for  any  particular  observation  point,  we  have  specified  a  technique  for 
deriving  the  optimum  frequency,  i.e.,  the  frequency  for  which  the  ampli¬ 
tude  of  the  electric  field  is  a  maximum. 

(2)  We  have  extended  the  work  of  Wait  [1952]  on  the  electric  fields 
produced  by  a  linear  current  source  of  finite  length  by  deriving  an  ex¬ 
pression  for  the  single  magnetic  field  component  (the  azimuthal  component) 
produced  by  such  a  source.  We  have  also  used  a  numerical  integration 
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technique  to  obtain  numerical  values  for  the  amplitudes  of  the  electric 
and  magnetic  fields  produced  in  sea  water  by  representative  finite 


sources. 

(3)  By  extending  one  end  of  the  finite  source  in  (2)  to  infinity, 
we  have  derived  new  expressions  for  the  electric  and  magnetic  fields 
produced  by  a  linear  current  source  of  semi- inf inite  length.  Once  again 
using  numerical  integration,  representative  numerical  values  for  the  two 
components  of  electric  field  have  been  c  -ained  around  the  end  (a  "point 
electrode")  of  this  semi-infinite  source  in  sea  water. 

(4)  By  combining  two  of  the  semi-infinite  sources  described  in  (3), 
we  have  derived  expressions  for  the  electric  and  magnetic  fields  produced 
by  a  linear  current  source  of  infinite  length  containing  a  gap  of  finite 
size.  Two  versions  of  this  source  are  considered:  (a)currents  of  equal 
ampl itudes  flowing  in  the  same  directions  in  both  semi -inf inite  segments, 
and  (b)  currents  of  equal  amplitudes  flowing  in  opposite  directions  in 
the  two  semi-infinite  segments.  Once  again,  representative  numerical 
values  have  been  computed  for  the  amplitudes  of  the  electric  fields  and 
the  fields  have  been  compared  with  those  produced  by  an  equivalent 
continuous  source  of  infinite  length. 

(5)  Finally,  we  have  reduced  the  size  of  the  finite  gap  in  (4) 
until  it  is  infinitesimally  small  and  shown  that  the  fields  produced  in 
the  vicinity  of  the  elementary  gap  are  equivalent  to  the  fields  produced 
by  a  continuous  linear  current  source  of  infinite  length  less  the  fields 
produced  by  an  elementary  current  source  at  the  location  of  the  gap. 

In  general,  we  have  made  as  few  assumptions  as  possible  in  deriving 
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the  above  expressions  and  numerical  data.  The  assumptions  that  are  made 
can  be  summarized  as  follows.  First,  the  displacement  current  is  neg¬ 
lected.  This  is  a  very  good  assumption  at  all  frequencies  up  to  about 
100  MHz  if  the  surrounding  medium  is  sea  water.  Second,  we  assume  that 
the  amplitude  of  the  alternating  current  is  the  same  throughout  the 
various  sources.  At  the  ULF/ELF  frequencies  that  are  considered,  this 
is  equivalent  to  assuming  that  the  current  sources  are  insulated  every¬ 
where  except  at  their  ends.  Third,  we  assume  that  the  shortest  distance 
between  any  observation  point  and  the  surface  of  the  sources  (which  are 
taken  to  be  circular  in  cross-section)  is  always  much  greater  than  the 
source  radius.  This  assumption  is  discussed  in  greater  detail  in  Appen¬ 
dix  A.  Finally,  as  we  have  pointed  out  many  times,  the  conducting  medium 
is  always  assumed  to  be  infinite  in  extent. 

VIII. 2  Discussion 

Our  first  objective  in  deriving  the  field  expressions  in  Chapters 
III,  V,  and  VI  for  sources  of  finite  and  semi-infinite  length,  and  of 
infinite  length  with  a  gap  of  finite  size,  was  to  investigate  the  possi¬ 
ble  desirable  characteristics  of  the  fields  near  their  open  ends  (i.e., 
point  electrodes)  for  undersea  communication  at  ULF/ELF  frequencies. 
Assuming  that  this  first  part  of  cur  study  gave  positive  results,  our 
second  objective  was  to  derive  enough  theory  and  numerical  data  to  pro¬ 
vide  the  basis  for  the  preliminary  design  of  undersea  communication 
arrays  that  would  have  the  capability  of  producing  measurable  ULF/ELF 
signals  throughout  limited  regions  of  the  ocean  without  the  cables  pro- 
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due i ng  the  signals  necessarily  being  located  in  the  regions.  To  a  large 
extent  our  objectives  have  been  achieved:  comparatively  large  electric 
and, to  a  lesser  extent,  magnetic  fields  can  be  produced  at  ULF/ELF  fre¬ 
quencies  in  the  vicinity  of  the  open  ends  of  the  sources  considered  in 
this  work,  and  measurable  ULF/ELF  fields  can  be  produced  in  and  surround¬ 
ing  the  region  of  a  moderately  large  gap  in  an  otherwise  continuous 
current-carrying  cable.  However,  it  is  important  for  us  to  point  out 
that  the  same  frequency  versus  range  constraint  applies  in  the  case  of 
the  electromagnetic  fields  produced  by  undersea  cables,  with  and  without 
gaps,  as  is  encountered  in  other  methods  of  undersea  communication  in¬ 
volving  electromagnetic  signal  propagation  in  sea  water:  large  ranges 
can  be  achieved  by  using  low  frequencies,  particularly  low  ULF  frequen¬ 
cies,  but  because  the  rate  of  data  transfer  is  directly  proportional  to 
the  signal  frequency, the  rate  of  data  transfer  over  the  large  ranges  may 
be  too  small  to  be  useful.  As  we  will  now  show,  this  frequency  versus 
range  constraint  effectively  limits  the  use  of  discontinuous  cable 
arrays  to  short  range  communication. 

Comparing  the  ULF/ELF  fields  produced  in  sea  water  by  cables  of 
finite  and  semi-infinite  length,  and  of  infinite  length  with  a  gap,  with 
the  fields  produced  by  a  long  continuous  current-carrying  cable  (our 
linear  current  source  of  infinite  length),  our  data  show  that  at  short 
distances  from  their  open  ends  i.e.,  distances  less  than  about  1  km,  the 
discontinuous  sources  can  produce  greatly  enhanced  electric  fields.  We 
have  not  carried  out  a  complete  computation  of  all  the  field  amplitudes 
that  can  be  produced  under  all  possible  circumstances,  but  the  available 
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data  indicate  that  the  increase  is  confined  largely  to  frequencies  in 
the  UIF  range  and,  within  this  range,  it  becomes  larger  as  the  frequency 
is  decreased.  At  0.001  Hz,  for  example,  the  increase  of  nearly  2000 
times  that  is  observed  at  a  perpendicular  distance  of  10  m  from  the  cen¬ 
ter  of  a  100  m  gap  in  an  otherwise  continuous  long  current-carrying  ca¬ 
ble  (Figure  VI. 2)  is  nearly  10  times  larger  than  the  increase  at  0.01  Hz. 
In  all  the  cases  we  have  examined  the  increase  begins  to  decline  rapialy 
with  distance  once  some  limiting  distance  has  been  exceeded.  For  the 
.wo  gap  systems  discussed  in  Chapter  VI  the  limiting  distance  is  of  the 
same  order  as  the  width  of  the  gap  (100  m) ,  a  result  that  may  be  indica¬ 
tive  of  a  general  relation  at  the  low  UIF  frequencies  under  considera¬ 
tion.  It  will  be  observed  that  the  increase  under  discussion  appears  to 
be  confined  to  the  electric  field,  as  may  be  seen  by  comparing  the  numeri 
cal  data  plotted  in  Figures  III. 9  and  III. 10  for  the  electric  and  magne¬ 
tic  fields  produced  by  a  cable  of  finite  length.  It  is  obvious  from 
these  and  the  other  results  we  have  derived  that  discontinuous  sources 
have  at  best  no  advantage  over  the  continuous  linear  source  fo>"  field 
generation  at  distances  much  above  1  km  and,  at  worst,  they  may  produce 
much  smaller  fields  than  the  continuous  source  for  some  geometries  and 
frequencies.  We  conclude  therefore  that  arrays  of  finite  and/or  semi¬ 
infinite  cables,  and  continuous  cable  with  gaps,  can  be  used  most  advan¬ 
tageously  for  undersea  communication  only  over  short  distances. 

The  possible  situations  where  such  communication  could  be  u.ed  are 
not  necessarily  trivial.  For  example,  the  ULF/EIF  electric  fields  around 
the  gap  between  two  point  electrodes,  or  near  a  single  point  electrode, 
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could  provide  a  low-data-rate,  non-acoustic  means  of  communication  with 
an  undersea  or  sea-floor  receiver  located  beneath  the  hull  of  an  elec¬ 
trode  or  cable/electrode  equipped  ship.  At  the  short  distances  under 
consideration  (distances  <  1  km)  there  is  likely  to  be  no  difficulty 
providing  communication  with  divers  at  frequen  ' in  the  lower  part  of 
the  voice  frequency  range  (50  -  3000  Hz);  indeeu,  a  commercial  unit  has 
been  designed  for  an  application  of  this  kind  [MacLeod,  1977].  As  a 
final  note  on  the  possn  c  comnunication  applications  of  the  discontin¬ 
uous  sources,  we  observe  that  long  range  is  not  always  a  desirable  fea¬ 
ture  cf  a  communication  system,  and  under  some  circumstances  the  limited 
ranges  for  signal  propagation  that  we  have  discussed  here  could  be  a 
decided  advantage. 

Turning  now  to  the  electromagnetic  fields  produced  at  large  dis¬ 
tances  from  the  sources  (distances  of  the  order  of  10  km,  or  greater), 
our  data  show,  as  we  have  commented  afo.e,  that  no  discontinuous  source 
produces  larger  fields  man  the  long  continuous  current-carrying  cable, 
and  most  of  the  sources  produce  smaller  fields.  Thus  long  continuous 
current-carrying  cables  appear  to  be  the  best  choice  if  long  distance 
undersea  communication  using  the  electromagnetic  fields  from  cable  sources 
is  the  goal.  From  the  results  given  in  Chapter  IV,  it  can  be  seen  that 
the  maximum  range  that  can  ne  achieved  with  an  alternating  current  of 
1000  A  with  present  electric  and  magnetic  ♦  i •  1  d  sensors  is  of  the  order 
of  0.1  -  4  km  for  frequencies  in  the  range  1H00  -  1  Hz  and  of  the  order 
o*  4  -  100  km  for  free  •■•"C  ier  in  *hf-  range  1  -  ''.r'r'  * <■_  assume  that 

the  cables  are  in  a  sea  of  ;  -te  extent  i>  ■  Its;  in 
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practice  they  will  have  to  be  located  on  the  sea  floor  and  sea-floor  and 
sea-surface  effects  will  influence  the  maximum  ranges  that,  ^an  be  achiev¬ 
ed.  but  we  would  not  expect  the  ranges  to  be  extraordinarily  different 
from  those  we  have  quoted  unless  the  sea  is  electrically  shallow  and  a 
significant  up-over-and-down  propagation  mode  exists  [Fraser- Smith  and 
Bubenik,  1979].  The  maximum  range  that  can  be  achieved  at  any  frequency 
can  be  made  larger  either  by  increasing  the  amplitude  of  the  current  or 
by  increasing  the  sensitivity  of  the  sensors  used  to  measure  the  electro¬ 
magnetic  fields.  However,  (1)  our  assumed  1000  A  current  is  already  sub¬ 
stantial  and  a  significant  increase  would  entail  large  increases  in  power, 
and  (?)  the  near-vertical  slopes  of  most  of  the  curves  in  Figures  IV. 1 
and  IV. 2  as  they  approach  the  lower  hori2ontal  axes  implies  that  even 
major  increases  in  sensor  sensitivity  would  only  have  a  small  effect  on 
maximum  range  (even  when  no  allowance  is  made  for  the  presence  of  a  back¬ 
ground  of  natural  noise).  The  range  can  always  be  increased  by  reducing 
the  ft equency,  but  the  lowest  frequencies  considered  in  Figures  IV. 1  and 
IV. 2  are  already  extraordinarily  low  by  present  undersea  communication 
standards,  which  severely  restricts  the  rate  of  data  transfer,  and  there 
is  the  further  disadvantage  that  the  background  noise  increases  'apidly 
with  decreasing  frequency  throughout  the  ULF  range.  We  conclude  that 
maximum  ranges  of  1  -  1C  km  are  likely  to  be  the  limit  for  a  single  long 
continuous  cable  and  that  undersea  communication  by  means  of  the  electro¬ 
magnetic  fields  from  UIF/EIF  currents  in  undersea  cables  is  necessarily 
restricted  to  special  applications  cf  limited  range  using  either  single 
long  cables  (transverse  range  1  -  10  km)  or  arrays  of  long  cables  (traes- 
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verse  ranges  of  a  few  tens  to  a  few  hundreds  of  kilometers,  depending  on 
the  number  of  cables). 


In  Appendix  0  we  consider  two  basic  arrays  of  long  cables:  a  series 
array,  and  a  parallel  array.  Assuming  once  again  that  our  infinite-sea 
ranges  are  comparable  to  those  that  would  be  computed  for  a  real  sea  with 
both  a  sea  surface  and  a  sea  floor,  the  results  derived  in  the  appendix 
indicate  that  a  series  array  of  10  cables  spaced  7.5  km  apart  and  driven 
by  a  1000  A  current  at  1  Hz  could  produce  a  measurable  magnetic  field  at 
a  distance  of  1  km  above  the  array  (observation  point  in  the  sea)  over  a 
transverse  distance  of  about  75  km  and  a  longitudinal  distance  limited 
essentially  to  the  length  of  the  cables.  The  natural  application  of  such 
arrays  is  to  provide  essential  communication  with  undersea  receivers 
throughout  limited  but  important  regions  of  the  sea. 

VIII. 1  Suggestions  for  Further  Work 

An  obvious  extension  of  our  work  would  be  to  include  a  sea  floor 
and  then,  in  addition,  a  sea  surface  in  our  theoretical  formulation. 

Such  an  extension  would  involve  a  major  new  theoretical  effort.  A  less 
obvious  extension  would  be  to  consider  the  use  of  pulsed  currents  in  the 
submerged  cables,  instead  of  the  continuous  alternating  currents  assumed 
in  this  work . 

The  above  extensions  are  important,  but  perhaps  the  most  important 
extensions  that  could  be  undertaken  is  a  feasibility  study  to  determine 
the  practical  utility  of  the  proposed  methods  of  undersea  communication. 
We  have  indicated  several  different  possible  methods  of  using  undersea 
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cables  for  communication,  but  we  have  necessarily  omitted  full  considera¬ 
tion  of  their  feasibility.  Instead  we  have  concentrated  on  making  avail¬ 
able  the  theoretical  expressions  and  some  of  the  numerical  data  that 
would  be  required  in  a  feasibility  study.  We  conclude  by  reiterating 
that  most  of  the  theoretical  expressions  and  numerical  data  we  have  de¬ 
rived  are  new  and  as  such  they  could  well  have  unexpected  implications 
for  undersea  communication.  It  would  certainly  be  a  departure  from  pre¬ 
sent  practice  to  communicate  up  to  a  submerged  receiver,  instead  of  down 
to  the  receiver,  and  to  have  a  signal-to-noise  ratio  that  increases  with 
depth  in  the  sea.  It  is  for  these  latter  reasons,  considered  in  relation 
to  the  growing  importance  of  reliable  undersea  communication  to  the  U.S. 
Navy,  that  the  suggested  feasibility  study  appears  particularly  appro¬ 
priate  at  this  time. 


77 


IX.  REFERENCES 


a 


Butterworth,  S.,  "The  distribution  of  the  magnetic  field  and  return 
current  round  a  submarine  cable  carrying  alternating  current.  - 
Part  2,"  Phil.  Trans.  Roy.  Soc.  London,  Ser.  A,  224,  141,  1924. 

Drysdale,C. V. ,  "The  distribution  of  the  magnetic  field  and  return  current 
round  a  submarine  cable  carrying  alternating  current.  -  Part  1," 
Phil,  Trans.  Roy.  Soc,  London,  Ser.  A,  224,  95,  1924. 

Dwight,  H.B.,  Tables  of  Integrals  and  Other  Mathematical  Data,  Fourth 
Edition,  336  pp.,  MacMillan,  1961. 

Fraser-Smith,  A.C.,  and  D.M.  Bubenik,  "ULF/ELF  electromagnetic  fields 
generated  above  a  sea  of  finite  depth  by  a  submerged  vertically* 
directed  harmonic  magnetic  dipole,"  Radio  Science,  14,  59,  1979. 

Kraichman,  M.B.,  "Handbook  of  Electromagnetic  Propagation  in  Conducting 
Media, "  Second  Printing,  U.S.  Government  Printing  Office,  Washing¬ 
ton,  D.C.,  1976. 

Lebedev,  N.N.,  Special  Functions  and  Their  Applications.  308  pp.,  Dover, 
N.Y.,  1972. 

Lowell,  H.H.,  Tables  of  the  Bessel-Kelvin  functions  ber,  bei,  ker,  kei, 
and  their  derivatives  for  the  argument  range  0(0.01)107.50,  Tech. 
Rept.  R-32,  NASA,  Lewis  Res.  Center.  Cleveland,  Ohio,  1959. 

MacLeod,  N.,  "Electric  diver  communication;  Non-acoustic  system  can 
operate  in  noisy  environments,"  Sea  Technology.  18,  21,  May  1977. 


79 


Panofsky,  W.K.H.,  and  M.  Phillips,  Classical  Electricity  and  Magnetism. 


494  pp.,  Addison-Wesley,  Mass.,  1962. 

Scarborough,  J.B.,  Numerical  Mathematical  Analysis,  Sixth  Edition, 
pp.  600,  Johns  Hopkins  Press,  Baltimore,  1966. 

Staff,  Computation  Lab.  of  Harvard  Univ.,  Tables  of  the  Generalized 
Exponential-Integral  Functions,  416  pp.,  Harvard  Univ.  Press, 
Cambridge,  Mass.,  1949. 

Sunde,  E.O.,  Earth  Conduction  Effects  in  Transmission  Systems, 

373  pp..  Van  Nostrand,  N.Y.,  1949. 

Von  Aulock,  W.,  "Propagation  of  electromagnetic  fields  in  sea  water," 
Tech.  Memo.  No.  140,  Dept,  of  the  Navy,  Bureau  of  Ships, 
Washington,  O.C.,  May  1948. 

Von  Aulock,  W. ,  "The  electromagnetic  field  of  an  infinite  cable  in  sea 
water,"  Tech.  Rept.  No.  106,  Dept,  of  the  Navy,  Bureau  of  Ships, 
Washington,  D.C.,  1  Sept.  1953. 

Wait,  J.R.,  "Electromagnetic  fields  of  current-carrying  wires  in  a 
conducting  medium,"  Canadian  J.  Phys. ,  30,  512,  1952. 

Wait,  J.R.,  "Radiation  from  current  distributions,"  pp.  167-200  in 
Electromagnetic  Radiation  from  Cylindrical  Structures,  Pergamon, 
N.Y.,  1959. 

Wait,  J.R.,  "Propagation  of  electromagnetic  pulses  in  a  homogeneous  con 
ducting  earth,"  Appl ied  Sci .  Res. ,  Sect.  B,  8,  213,  1960. 


80 


Wright,  C., 


Electromagnetic  Phenomena,  Ed.  G.W.  Wood,  216  pp. ,  December  1953. 

Young,  A.,  and  A.  Kirk,  Bessel  Functions,  Part  IV;  Kelvin  Functions, 
Royal  Soc.  Math.  Tables,  Vol.  10,  pp.  97,  Cambridge  Univ.  Press, 


APPENDIX  A 


THE  INFINITELY-THIN  LINEAR  CURRENT  SOURCE  APPROXIMATION 


In  Section  II  we  found  that  the  inhomogeneous  wave  equation  (Equa¬ 
tion  1 1. 9)  had  a  Hert2  vector  solution  of  the  form 
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in  an  unbounded  homogeneous  conducting  medium.  Here  r  is  the  distance 
between  the  volume  element  dv  and  the  observation  point  P.  Note  that  we 
use  primed  cylindrical  coordinates  (o',  b‘,  z')  for  source  point  and  un¬ 
primed  coordinates  (o,  z)  for  observation  points.  Our  object  in  this 
Appendix  is  to  derive  a  version  of  the  above  Hertz  vector  solution  for  a 
linear  current  source.  As  we  will  show,  it  is  possible  to  obtain  an 
approximate  solution  in  which  the  volume  integral  is  replaced  by  a  line 
integral  along  the  axis  of  the  source. 

Figure  A.l  shows  the  cylindrical  linear  current  source  oriented  in 
the  z-direction  and  extending  from  z‘  =  ^  to  z'  =  Its  radius  is  pc, 
and  it  carries  current  of  density  Jzexp(iujt).  From  Equation  (A.l),  the 
Hertz  vector  for  this  case  has  only  a  z-component  which  can  be  written  as 
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We  will  assume  that  the  current  density  is  a  function  of  z'  only. 
We  will  then  consider  the  following  two  conditions: 

(i)  p  >>  p‘  for  all  p',  so  that  (p  +  o')  ^  p. 

(if)  (z-z1)2  »  (p+p')^,  so  that  (p  +  p')2  +  ( z-z  * ) 2  =  (z-z’)c  . 

These  conditions  have  the  following  physical  significance.  First, 
condition  (i)  implies  that  the  p-component  of  the  observation  point  must 
be  much  larger  than  the  radius  of  the  cylindrical  linear  current  source. 
Second,  condition  (ii)  implies  that  the  shortest  vertical  distance  from 
the  observation  point  to  the  cylindrical  linear  current  source  is  much 
larger  than  both  the  p-component  of  the  observation  point  and  the  radius 
of  the  source.  Note  that  the  two  conditions  do  not  necessarily  apply 
simul taneously. 

If  the  first  condition  is  true,  it  follows  that  r  s  [p  +  (z-z‘)  ]  . 
If  the  second  condition  is  true,  it  follows  that  r  =  jz-z'j.  If  any  of 
these  conditions  is  true,  the  distance  r  is  no  longer  a  function  of  p’ 
and  o'  and  the  Hertz  vector  can  be  rewritten  as 
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where  the  two  integrals  on  the  right  hand  side  give  irpc  .  Since 
J z ( z  ‘ )  =  I z ( z ’ ) ,  it  follows  that 
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Figure  A.1.  Cylindrical  antenna  of  finite  length  (f?  -  l] )  and  finite 
cross-section  (radius  r>c )  carrying  a  current  with  a  uniform  density 


A  small  volume  element  dv  of  the  antenna  is  also  shown  in  the 
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In  our  derivations,  we  have  always  used  Equation  (A. 5)  for  the  Hertz 
vector.  This  means  that  all  our  field  expressions  are  valid  as  long  as 
one  of  the  above  conditions  (i)  and  (ii)  are  satisfied,  i.e.,  as  long  as 
the  shortest  distance  from  the  observation  point  to  the  cylindrical 
linear  current  source  is  much  larger  than  the  radius  of  the  cylindrical 
linear  current  source  itself.  In  other  words,  we  assume  that  the  linear 
current  source  is  infinitely  thin. 


append: x  b 


INTRODUCTION  OF  THE  KELVIN  FUNCTION, 


In  Section  IV,  Equation  (IV. 2),  we  have  an  integral  of  the  form 


•0  +  i)u 


(B.l ) 


2  2  ^ 

where  u  =  (a  +  x  ) ,  which  we  want  to  express  in  terms  of  known 


functions. 


This  integral  can  also  be  written  as 


30  -(l+i)u 


(B.2) 


Defining  u  «  at,  this  integral  can  be  rewritten  as 


;i  =x 


30  e-(l  +  i)at 

HT7 

V  1 


(B.3) 


From  Lebedev  [1972;  see  p.  119],  we  have  the  following  integral 
expression  for  the  n-th  order  modified  Bessel  function  of  the  second 
kind: 


Kn(z)  =  fGk)  ^  f,  e"Zt  (B-4) 


provided  Re  z  >  o  and  Re  n  >  -4.  Here,  T  is  the  Ga.nma  Function  defined 


r  (n)  =  I  yn 
'o 


e'^dy. 


( B.  5) 


Choosing  n  =  0  and  knowing  that  F  (^)  =  <Jr,  ,  Equation  (6.4)  becomes 


r 


KQ(z)  =  /  -=■-  dt. 


(B.6) 


‘1  W  t2-] 


Substituting  z  =  a(l+i)  in  Equation  (B.6)  and  comparing  it  with  Equation 
(B.3) ,  we  obtain 


Ko  (  a<1  +  i))  =  '}■  (5.7) 

In  other  words,  the  integral  I-|  is  a  zero-th  order  modified  Bessel  func¬ 
tion  of  the  second  kind  of  complex  argument.  Knowing  this  fact,  we 
can  now  proceed  to  write  I-j  in  terms  of  known  functions. 

Since  y  =  (l+i)3  and  a  =  Be,  we  have 


K0  (  *(1  +  i)  )  =  K0M 


and  from  these  relations  it  follows  that 


-1 


■  V») 


/ 


u 


dx . 


We  also  know  from  Dwight  [1961]  that 


( B  .8) 


(B.9) 
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i  0  Kn  |  a  JT  )  -  kern  a  +  i  kein  a  ,  (B.10) 


where  ker  a  and  kei  a  are  Kelvin  functions  of  the  second  k  .d  of  order 
n  n 

n  [Young  and  Kirk,  1964].  The  left  hand  side  of  Equation  (B.10)  can 
also  be  written  as 

Kn(a  '/7’)=  (6'n> 

Letting  n  =  0  and  following  equations  (B.8)  -  (B.ll),  we  have 

I,  =  Ko(yp)  =  kero(jTa)  +  i  kei0(/Fa)  .  (B.12) 


We  have  now  expressed  the  integral  I j  in  terms  of  known  Kelvin 
functions.  If  we  decompose  in  equation  (B.9)  into  its  real  and 
imaginary  parts,  we  obtain 


=  K  (yp)  « 

o'  ■  ‘ 


/ 


e  cos  u 


dx  - 


oo 

■  / 


-U. 


sin  u 


u 


dx  ,  (B.l 3) 


then,  comparing  equations  (B.12)  and  (B.13),  it  follows  that 


(B.  14) 
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If  we  let  the  upper  limits  of  the  integrals  given  by  Equation  (III. 8) 


go  to  infinity  and  compare  the  results  with  the  integrals  in  Equation 
(B.14) ,  we  can  write 


kero  (/2~a)  -  f  7Tdx  *  Ec  “3  • 

0 

(B. 1 5) 

ke1o  =  -  Es  Ca»  “]  - 


Similarly,  if  we  also  let  the  upper  limits  of  the  integrals  given 
by  Equation  (III. 16)  go  to  infinity  and  compare  the  results  with  the 
derivatives  of  the  Kelvin  functions  given  by  Equation  (B.14),  we  obtain 


kero  a 3  =  (  McCa,  °°]  +  Ms[a,  ®]  +  Nc [a ,  «]  j 

y  /  n 


(B.16) 


ke’o  a^  s  ~nT  (  McC*’  "  Ns^a'  ) 


Since  we  also  know  that 


ker:  <JTa>  *  -jy  (kar0  </2  a)  -  k*<i  (J2  a>)  • 


(B. 1 7) 


ke1i  ('^”a)  '  jy  (keri  (/^a)  *  ke1i  • 

we  can  substitute  Equations  (B.16)  into  Equations  (3.17)  and  obtain 


ker1  (JT  a) 
kei'i  (JTa) 


~T  (  2Ms  °°]  4  Nc  [a ,  oo]  +  N$  [a,  oo]  ^  , 

( B. 18) 

-f-  2Mc  [a,  »]  +  Nc  [a,  «]  -  N$  [a,  °°]j  . 
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APPENDIX  C 


FIELDS  PRODUCED  BY  A  DC  CURRENT 


In  this  appendix  we  summarize  the  electric  and  magnetic  field 
expressions  for  dc  currents  flowing  in  linear  current  sources  of 
infinitesimal,  finite,  semi-infinite,  and  infinite  length. 

1 .  Current  Element  Source 

The  fields  produced  by  linear  current  element  of  length  d*.  and 
carrying  a  dc  current  I  are  given  by 


E  « 


3Idt  sine  cose 
4rar3 


E7  =  (3  cos26  -  1 ) 

4TOr 


H 


IdL  sine 


*  47ir2 


and 


(C.l) 


2 .  Linear  Current  Source  of  Finite  Length 

The  fields  due  to  a  linear  current  source  of  finite  length  carrying 
a  dc  current  I  and  extending  from  to  are 
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E  =  -fii- 
P  4tto 


1  _  J_ 

TP  ry 


E  =  — 

tZ  4'n'O 


Z-l2  z~l-\ 


r  3  r  3 

r2  r1 


(C.2) 


and 


H  = 

4>  4ti  p 


z-ly  z-^2 


with  E  =  H  =  H  =  0  . 

<f  P  z 


Now  consider  a  perpendicular  line  passing  through  the  center  of  the 


source.  Along  this  line,  E  =  0  and  we  have 

P 


E  ^  _lIL_ 

3  ’ 

47,orJ 


and 


H„  = 


IL 


4npr  ’ 


where  L  =  -  Z-j  is  the  length  of  the  source. 

If  we  normalize  all  the  distances  to  the  length  L  of  the  source, 
we  obtain 


L2E  =  — 1 1— 


4nori 


and 


(C.  3) 


LH.  = 


t  4^rpLrL 


where  p,  =  p/L  and  r.  =  r/L. 


In  Figure  C.l,  we  have  plotted  the  variation  of  the  magnitude  of 

p 

the  normalized  E2  field  (i.e.,  EzL  )  with  normalized  distance  p^. 


3.  Linear  Current  Source  of  Infinite  Length 
Given  a  linear  current  source  of  infinite  length  carrying  a  dc 
current  I,  all  the  field  components  are  zero  except  for  the  magnetic 

field  component  H.,  where 

$ 


H  = 

<p  2ttc 


4.  Linear  Current  Source  of  Semi-Infinite  Length 
In  the  case  of  a  semi-infinite  linear  current  source  extending  from 
-«>  to  and  carrying  a  dc  current  I,  we  have 


3  * 


z  4-or23  2 


(C.5) 


I  ,  z-£; 

H4>  "  4^p  '  r? 


E.  =  H  =  H  =  0  . 
$  P  z 


5 .  Linear  Current  Source  of  Infinite  Length  with  a  Finite  Gaf 
Assuming  the  direction  of  current  f  1  c"-'  is  the  same  in  the  two 
segments  of  the  source  (Figure  VI. 1(a)),  we  have 
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(C.6) 


and 


Again, E  =H  =  H  =0. 

i  P  2 

These  equations  are  the  negative  of  Equations  (C.2),  which  means 
that  in  the  dc  case,  the  gap  system  produces  fields  of  the  same  magni¬ 
tude  but  pointing  in  opposite  directions  compared  with  a  linear  current 
source  of  finite  length  carrying  the  same  current  I  as  flows  through 
the  gap  and  occupying  the  same  position  as  the  gap.  Thus  if  the  gap 
length  is  L  (i.e.,  =  L),  the  plot  of  the  magnitude  of  the  electric 

field  along  the  perpendicular  line  passing  through  the  center  of  the  gap 
will  be  exactly  the  same  as  the  one  shown  for  a  linear  current  source  of 
finite  length  in  Figure  C.l. 

If  the  currents  in  the  two  semi-infinite  segments  are  directed 
towards  each  other  (See  Figure  VI. 1(b)),  instead  of  in  the  same  direction 
as  above,  the  electric  and  magnetic  field  expressions  become 
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Figure  C.l.  Variation  with  normalized  perpendicul  ar  distance  p/L  of  the 
total  electrical  fields  produced  in  sea  water  by  linear  current  sources 
carrying  dc  currents  of  1000  A.  One  curve  (circles)  applies  to  both  a 
linear  current  source  of  finite  length  L  and  to  a  linear  current  source 
of  infinite  length  having  a  finite  gap  of  length  L  with  the  current 
flowing  in  the  same  direction  in  each  semi-infinite  segment.  The  other 
curve  (squares)  corresponds  to  a  linear  current  source  of  infinite 
length  having  a  finite  gap  of  length  L  with  the  current  flowing  in 
opposite  directions  in  each  semi -inf ini te  segment.  The  perpendicular 
axis  along  which  the  distance  variations  are  shown  passes  through  the 
central  point  of  each  source. 
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and 


H  =  ~- 

<J>  4trp 


2-S.i 


Z-2.„ 


If  we  now  restrict  ourselves  to  the  fields  produced  along  the 
perpendicular  line  passing  through  the  center  of  the  gap  of  finite 
length  L  and  normalize  all  the  distances  with  L,  we  obtain 


2ncr[_' 


E 

z 


=  0  , 


(C.8) 


H*  =  0  . 

T 

2 

The  variation  of  the  magnitude  of  the  normalized  Ep  field  (i.e.,  EpL  ) 
with  normalized  distance  for  this  case  is  plotted  in  Figure  C.l. 

If  we  now  compare  the  electric  field  variations  shown  in 
Figure  C.l,  we  find  that  the  electric  field  produced  by  the  gap  system 
with  opposite  current  flow  shown  in  Figure  VI. 1(b)  declines  less  rapidly 
with  distance  (once  r  exceeds  L/2)  than  the  electric  field  produced 
eitner  by  the  gap  system  with  currents  flowing  in  the  same  direction  in 
each  semi -infinite  segment,  as  shown  in  Figure  VI. 1(a),  or  by  the  linear 
current  source  of  length  l.  The  respective  slopes  of  the  two  curves  are 
-20  dB/decade  and  -30  dB/decade  at  large  distance. 
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APPENDIX  D 


FIELDS  PRODUCED  BY  ARRAYS  OF  LONG  CABLES  SUBMERGED 
IN  A  CONDUCTING  MEDIUM  OF  INFINITE  EXTENT 

As  discussed  in  Chapter  VIII,  the  data  displayed  in  the  figures  in 
Chapter  IV  suggest  that  a  single  submerged  cable  cannot  produce  electric 
or  magnetic  fields  with  the  right  frequency/amplitude  characteristics 
for  successful  undersea  communication  once  the  range  exceeds  some  value 
that  will  typically  lie  in  the  range  1  -  10  km  (the  actual  value  will 
depend  on  a  number  of  practical  details;  some  of  the  most  important  of 
these  details  are  the  frequency,  the  amplitude  of  the  current  in  the 
cable,  the  amplitude  of  the  signal,  the  range,  the  background  noise 
level,  and  the  sensitivity  of  the  receivers).  However,  there  appears  to 
he  no  reason  why  the  range  cannot  be  extended  more  or  less  indefinitely, 
in  principle,  by  the  use  of  arrays  of  long  cables.  The  problem  of  range 
is  then  converted  to  a  problem  of  cost  and  of  power  availability.  In 
this  appendix  we  briefly  consider  some  of  the  requirements  for  these 
arrays  and  the  ranges  that  might  be  achieved. 

As  a  first  example  of  an  array  of  long  cables,  consider  the  config¬ 
uration  of  cables  (length  L)  shown  in  the  top  panel  of  Figure  D.l.  We 
will  refer  to  this  configuration  as  a  series  array,  and  we  will  assume 
that  L  is  so  large  that  we  can  treat  each  cable  in  the  array  as  a  linear 
current  source  of  infinite  extent  and  thus  compute  the  electric  and  mag¬ 
netic  fields  they  produce  from  the  expressions  and  data  given  in  Chapter 
IV.  We  also  assume  that  the  spacing  s,  although  small  compared  with  L, 
is  large  enough  for  each  cable  to  be  treated  as  being  independent  of  its 
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neighbors.  A  second  example  of  an  array  of  long  cables,  a  parallel 
array,  is  shown  in  the  bottom  panel  of  Figure  0.1;  we  will  make  the  same 
assumptions  for  this  array  as  those  we  have  just  detailed  for  the  series 
array.  We  will  suppose  that  the  origin  of  a  3-dimensional  cartesian 
coordinate  system  (x,y,z)  is  located  at  the  center  of  each  array  and 
that  the  planes  of  the  two  arrays  correspond  to  the  x,y  planes  of  the 
coordinate  systems,  with  the  y  axis  taken  parallel  to  the  cables  as 
shown  in  the  figure. 

Cross-sections  of  the  two  arrays,  taken  along  the  x-axis  in  Figure 
0.1,  are  shown  in  Figure  D.2.  Here  the  cables  are  designated  C  n,  C 
...  C_-j,  C+^ ,  ...  Cn  according  to  their  position  relative  to  the  origin 
of  the  coordinate  system.  Considering  the  series  array  first,  it  is 
easily  seen  that  the  total  electric  field  at  an  observation  point  loca¬ 
ted  at  a  height  h  above  the  origin,  i.e.,  at  the  point  (0,0, h)  is  always 
zero.  However,  the  magnetic  fields  produced  by  each  cable  add  up  vec- 
torially  and  give  a  total  magnetic  field  that  is  directed  along  the  ?.- 
axis.  The  first  two  vector  contributions  to  this  field,  the  magnetic 
fields  B+.|  and  B_-j  produced  by  the  two  nearest  cables,  are  shown  in  the 
figure.  The  electric  and  magnetic  fields  produced  by  the  parallel  array 
differ  substantially  from  those  produced  by  the  series  array.  Consider¬ 
ing  the  fields  produced  at  the  point  (0,0, h)  once  again,  we  find  that  the 
parallel  array  does  produce  an  electric  field,  unlike  the  series  array, 
and  the  field  is  directed  along  the  y  axis.  Further,  instead  of  being 
directed  along  the  z-axis,  i.e.,  in  a  direction  perpendicular  to  the 
plane  of  the  array,  the  magnetic  field  produced  by  the  parallel  array  is 
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Figure  b.l.  Two  different  arrays  of  parallel  long  cables  (length  =  L, 
spacing  =  s;  s«'L)  lying  in  the  same  plane  and  all  carrying  the  same 
current  I.  In  the  top  panel  (a)  we  show  a  series  array,  where  the 
current  in  a  particular  cable  flows  in  the  opposite  direction  to  the 
current  in  the  adjacent  cables.  In  the  bottom  panel  (b)  we  show  a 
parallel  array,  where  the  current  in  all  the  cables  flows  in  the 
same  direction.  Note  that  the  total  current  flowing  through  the 
series  array  is  I,  whereas  the  total  current  flowing  througn  the 
parallel  array  is  2nl,  where  2n  (n  =  1,2,3...)  is  the  number  of  cables 
in  each  array. 
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directed  along  the  x-axis,  i.e.,  in  a  direction  parallel  to  the  plane  of 
the  array.  These  various  properties  of  the  fields  do  not  apply  general¬ 
ly,  but  they  illustrate  the  large  differences  that  can  occur  in  the 
fields  produced  by  the  two  arrays. 

In  considering  arrays,  our  primary  interest  is  to  produce  measurable 
electric  and/or  magnetic  fields  throughout  as  large  a  volume  of  the  sea 
as  possible,  with  the  least  number  of  cables.  Cr,  since  the  sea  is  com¬ 
paratively  limited  in  depun  compared  with  horizontal  extent,  we  wish  to 
Droduce  measurable  fields  throughout  the  sea  over  as  large  an  area  as 
possible.  Clearly  the  cables  in  the  arrays  (which  we  now  assume  are  lo¬ 
cated  on  the  sea  floor,  ai trough  nis  is  not  a  necessary  condition  for 
the  use  of  the  arrays  we  have  been  considering)  must  be  spaced  far  apart, 
with  the  actual  spacing  being  determined  by  the  fields  produced  at  points 
midway  between  the  cables,  where  the  attenuation  of  the  fields  is  rela¬ 
tively  the  strongest  and  where  the  field  amplitudes  are  in  most  cases  a 
minimum  for  a  given  height  h  above  the  plane  of  the  arrays.  To  illus¬ 
trate  the  amplitudes  of  the  fields  that  can  be  produced,  we  have  taken 
the  series  array  shown  in  Figure  0.1(a)  and  D.?(a)  and  computed  the  vari¬ 
ation  of  the  amplitude  of  the  total  magnetic  field  with  height  h  along 

two  axes  perpendicular  to  the  plane  of  the  array:  the  z-axis,  which 

passes  through  a  point  (the  origin)  located  midway  between  two  cables, 

and  an  axis  passing  through  one  of  the  cables.  Because  our  interest  is 
in  t>  j  fields  produced  by  an  array  of  widely  spaced  cables,  we  assume 
that  the  distance  s  between  the  cables  is  30  skin  depths,  which  means 
that  we  can  reasonably  neglect  the  fields  produced  along  our  two  chosen 
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Figure  D.2.  Cross-sectional  views  of  the  two  arrays  of  long  cables 
shown  in  Figure  0.1.  At  top  is  the  series  array,  where  the  direction 
of  the  current  in  any  cable  is  opposite  to  the  direction  of  the  current 
in  adjacent  cables;  at  bottom  is  the  parallel  array,  where  the  direc¬ 
tion  of  the  current  is  the  same  in  all  cables.  The  magnetic  fields 
B_i  and  B+]  produced  at  the  observation  point  (o,o,h)  by  the  two  cables 
C-i  and  C+i  are  shown  for  each  array. 


axes  by  all  the  cables  in  the  array  except  for  the  one  or  two  cables 

nearest  to  the  axes.  Equation  (IV .6) is  used  to  obtain  the  magnetic  fields 

produced  by  these  cables. 

Figure  0.3  shows  the  results  of  the  computations.  The  data  are 
presented  in  parametric  form,  so  that  they  apply  at  any  ULF/ELF  frequency, 
and  it  is  immediately  apparent  by  comparing  the  two  curves  in  the  figure 
that  the  total  magnetic  field  produced  along  the  axis  passing  through 
one  of  the  cables  is  very  much  greater  than  the  field  produced  on  the 
axis  passing  through  a  point  midway  between  two  of  the  cables.  Using 
the  results  for  the  latter  axis  and  assuming  that  the  minimum  value  of 
the  magnetic  field  that  can  be  measured  is  10’^  pT,  we  find  that  for  a 
100  Hz,  1000  A  current  in  the  cables  the  magnetic  field  can  be  measured 

up  to  h  =  200  m  away  from  the  array,  and  for  1  Hz,  1000  A  current  in  the 

cables  the  magnetic  field  can  be  measured  up  to  h  -  1  km  away  from  the 
array. 

The  above  field  strength  considerations  are  perhaps  the  most  impor¬ 
tant  in  evaluating  the  possible  undersea  communication  performance  of 
cable  array.  However,  there  are  a  few  other  considerations  that  are 
also  important,  which  we  list  before  concluding  this  appendix.  First, 
because  communication  is  to  take  place  from  the  bottom  upwards  in  the 
sea,  it  may  not  be  necessary  for  measurable  electric  or  magnetic  fields 
to  be  produced  at  the  surface.  Second,  because  most  and  perhaps  all 
natural  ULF/ELF  noise  in  the  sea  comes  from  sources  above  the  surface, 
the  si gnal -to-noise  ratio  for  communications  via  undersea  cable  arrays 
should  increase  with  depth.  Further,  because  the  electric  and  magnetic 


01  234  56789  10  11  12 

h  /& 


Figure  0.3.  Variation  with  perpendicular  distance  h  (measured  in  skin 
depths  6)  of  the  amplitude  of  the  total  magnetic  field  produced  in  a 
conducting  medium  above  the  plane  of  the  series  array  shown  in  Figs. 
D.l(a)  and  0.2(a).  It  is  assumed  that  the  cables  in  the  array  carry 
a  1000  A  alternating  current  and  that  they  are  spaced  306  apart.  The 
top  curve  (circles;  left  hand  scale)  shows  the  magnetic  field  produced 
along  an  axis  perpendicular  to  the  plane  of  the  array  and  passing 
through  one  of  the  cables;  the  bottom  curve  (crossed  circles;  right 
hand  scale)  shows  the  magnetic  field  produced  along  an  axis  perpendi¬ 
cular  to  the  plane  of  the  array  and  passing  through  a  point  located 
midway  between  any  two  of  the  cables.  The  magnetic  field  direction 
for  the  top  curve  is  along  the  x-axis,  i.e.,  parallel  to  the  plane  of 
the  array  and  perpendicular  to  the  direction  of  the  cables,  and  for 
the  bottom  curve  the  magnetic  field  direction  is  along  the  z-axis, 
i.e.,  perpendicular  to  the  plane  of  the  array. 
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fields  from  these  natural  sources  are  predominantly  horizontally  directed 
in  the  sea,  it  may  be  possible  to  discriminate  against  the  natural  noise 
by  generating  and  measuring  vertically-directed  fields  in  the  sea. 
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